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Abstract 

' We present a comprehensive study of five-dimensional brane-world models for 

neutrino physics based on flat compactifications. Particular emphasis is put on 
^ . the inclusion of bulk mass terms. We derive a number of general results for such 

H ' brane-world models with bulk mass terms. In particular, in the limit of small 

brane-bulk couplings, the electroweak eigenstates are predominantly given as 
a superposition of three light states with non-trivial small admixtures of bulk 
states. As a consequence, neutrinos can undergo standard oscillations as well 
as oscillation into bulk Kaluza-Klein states. We use this structure to construct 
a specific model based on Z2 orbifolding and bulk Majorana masses which is 
compatible with all observed oscillation phenomena. The solar neutrino deficit 
is explained by oscillations into sterile bulk states while the atmospheric neu- 
trino deficit is due to ~ Ur oscillations with naturally maximal mixing. In 
addition, the model can accommodate the LSND result and a significant neu- 
trino dark matter component. We also analyze the constraints from supernova 
energy loss on neutrino brane-world theories and show that our specific model 
is consistent with these constraints. 



1 Introduction 



Two recent ideas, namely the brane-world idea |, ^, ^, ^ and the possibihty of 

having large gravitation-only additional dimensions |9|, | , |7 ] , may significantly change the 
relation between string- or M-theory and low-energy particle physics. Not only do these 
ideas motivate new directions in particle physics model building but they also provide 
new generic structures within string theory which may be experimentally testable. In 
this context, an important role is played by bulk particles coupling to standard model 
particles only gravitationally. Particularly relevant are the Kaluza-Klein modes of the 
higher-dimensional graviton Q and higher-dimensional bulk fermions. In this paper, we 
focus on the latter possibility of bulk fermions. They provide candidates for right-handed 
neutrinos and might, therefore, play an important role in neutrinos physics. 

The potential relevance of these particles for neutrino physics has first been first pointed 



out and analyzed in Ref. |Tn, 11, O]. Phenomenological analysis have been performed 



in |l^-|2g]. The importance of bulk masses has been pointed out in |20|, where an explicit 
6-dimensional example with Dirac bulk masses has been presented. The phenomenological 
implications of a five-dimensional brane-world model with bulk Dirac-masses for solar, 
atmospheric, short-baseline and supernova oscillations have been studied in p2| . 

The main goal of this paper is to present a comprehensive phenomenological study of 
five-dimensional brane-world models for neutrino mass and mixing. Particular emphasis 
will be put on the model-building options that arise from the various types of bulk mass 
terms. We will focus on compactifications on a flat metric as well as on bulk mass terms 
without any explicit dependence on the additional coordinate. Furthermore, we assume 
that all right-handed neutrinos propagate in the fifth dimension. Within this general 
setting we will present a number of exact results for the mass spectrum as well as the mixing 
angles. Furthermore, based on these general results, we will give an explicit example based 
on the orbifold 3^/^,2 and Majorana bulk-masses related to the S*^ model with Dirac 



masses studied in Ref. [22]. As we will see, this model is consistent with the oscillation 
data presently available. In particular, we show that the solar neutrino deficit can be 
explained by small-mixing angle oscillations of v^. into a tower of sterile neutrinos. Unlike 
such oscillations into a single sterile neutrino this option is not disfavored by the recent 
SuperKamiokande results. In addition, the model naturally leads to a maximal mixing 
angle in the - sector and can explain the atmospheric neutrino results. At the 
same time it contains two neutrinos in the eV range that can provide a cosmologically 
significant source of dark matter. It also has room to accommodate the LSND signal and 
in general predicts v^. <-> oscillations with Am^ > Am^rpj^ and a small amplitude. 
The alternative model-building option, namely large mixing angles between standard and 
sterile Kaluza-Klein neutrinos, is disfavored by the requirement of not having significant 



departures from Standard Model processes [|l^. In particular, this disfavors the possibility 



that sterile bulk states play a significant role in atmosperic neutrino oscillations. 

In the context of string- and M-theory, bulk fermions arise as superpartners of gravita- 
tional moduli, such as, for example, radii of internal spaces. Given this origin, the existence 
of bulk fermions is practically unavoidable in any supersymmetric string compactification 
and, therefore, represents a quite generic feature of string theory. This, in our opinion, 
constitutes the most likely origin of such particles within a fundamental theory and, at 
the same time, provides the best theoretical motivation to study brane-world neutrino 
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physics. In the context of four-dimensional effective actions from string theory, moduhni 
as candidates for right-handed neutrinos have been proposed in Ref. ||2^. The phenomeno- 
logical relevance of such bulk fermions for neutrino physics is basically controlled by three 
different features of the theory. Specifically, these are the masses of the bulk fermions, 
the radii of the additional dimensions and the size of the brane-bulk coupling between 
bulk fermions and standard model neutrinos. Let us discuss these three features and their 
relation to string- and M-theory separately. 

So far, brane-world model building in the context of neutrino physics has been mostly 
focusing on exactly massless bulk fermions (see, however, Ref. |l^ where bulk mass 
terms have been considered). However, bulk mass terms constitute an important model- 
building option which should be taken into consideration. In fact, unless forbidden by 
specific symmetries, bulk mass terms should not be ignored. Of course, bulk fermions 
can be candidates for sterile neutrinos only if these masses are sufficiently small Q. One 
should, therefore, worry about the origin of such small mass scales. Within perturbative 
string theory, moduli fields constitute flat directions and, consequently, have vanishing 
mass. The same is true for their fermionic superpartners as long as supersymmetry is 
unbroken. String theory, therefore, provides a generic reason why bulk fermion masses 
may be much smaller than the string scale. Eventually, non-perturbative effects have to 
be taken into account in order to break supersymmetry and stabilize the moduli fields. 
Those will generate masses for the moduli field and, most likely, masses for the bulk 
fermions as well. Hence, the inclusion of small bulk fermion masses is quite well motivated 
from the viewpoint of string theory. It is one of the central points of this paper, to take 
this insight seriously and include the most general structure of bulk mass terms into our 
analysis. 

The specific size of these mass terms as well as the size of the additional dimension 
depend on non-perturbative effects and are largely uncertain from a theoretical perspec- 
tive. Therefore, in the general part of the paper, we will not focus on any specific choice 
for these scales. Of course, the tower of bulk fermions will only be directly relevant for 
oscillations if the bulk masses and the scale of the additional dimension are sufficiently 
small. For the explicit example we will focus on such a case with small bulk masses and 
a large fifth dimension close to the experimental upper bound. 

Finally, we should address the possible origin of the brane-bulk mass terms mixing bulk 
fermions and standard neutrinos. In fact, such couplings may arise quite generically, for 
example, in the context of Horava-Witten theory Q. In the 11-dimensional formulation 
of this theory, the 10-dimensional Yang-Mills theories on the boundaries contain terms 
bilinear in the gravitino (giving rise to bulk fermions in lower dimensions) and the gauginos 
(giving rise to the standard model fermions). These are precisely the terms needed to 
generate the desired brane-bulk mixing. From the perspective of effective four-dimensional 
supergravity theories such terms can be understood as arising from higher-dimensional 
non-renormalizable operators involving moduli and standard model fields with appropriate 
VEVs inserted. It would certainly be interesting to analyze this in more detail. However, 
for the purpose of the present paper, we content ourselves with these general remarks. 

We would like to point out that there exist models, particularly in the context of string 



^ As we will see in Section 2, only the smallness of Lorentz-invariant mass terms needs to be accounted 

for. 



2 



compactifications, where, in addition to bulk right-handed neutrinos, one also has right- 
handed neutrinos on the brane. Such models are problematic for values of the fundamental 
scale far below the GUT scale. This is because a see-saw mechanism on the brane, with 
the right-handed neutrino masses being set by this low fundamental scale, will typically 
not suppress the associated contribution to the left-handed neutrino masses enough. In 
this paper, we will not consider models with such singlet neutrino states on the brane. 

In the next Section we present the general structure of the five-dimensional brane- world 
model along with the associated four-dimensional effective action. The Section concludes 
with a discussion of the various symmetries that can be imposed on this brane- world theory 
for model-building purposes. The spinor properties used in this Section are summarized in 
Appendix Section 3 contains a general analysis of the structure of masses and mixing 
and a discussion of the various resulting model building options. This discussion is mainly 
based on a perturbative diagonalization of the mass matrix assuming small brane-bulk 
mixing. The Section is complemented by Appendix ^ which contains some of the more 
technical aspects related to the exact diagonalization of the mass matrix. In Section 4 the 
constraints from energy loss in supernovae are discussed. Section 5 describes a specific 
model with a bulk Major ana mass which results in a small mixing between electroweak 
eigenstates and bulk states. It is shown that, in the context of this model, the small- angle 
MSW solution is a viable solution for the solar neutrino problem. At the same time, the 
atmospheric oscillation phenomena as well as the LSND signal can be understood. 

2 General framework 

In this Section, we will set up our formalism and present the general structure of 
five-dimensional neutrino brane-world models. Subsequently, we perform a Kaluza-Klein 
reduction to obtain the associated four-dimensional effective theory which will be the 
starting point for our analysis of neutrino masses and mixing. 

2.1 Five-dimensional brane-world models of neutrino physics 

We consider brane world theories with a five-dimensional bulk and coordinates denoted 
by x° where a, /3, ••• = 0,1,2,3,4. Spinors in this five-dimensional bulk theory provide 
candidates for right-handed neutrinos and are, therefore, of particular importance in our 
context. As mentioned above, within the framework of string- or M-theory compactifi- 
cations the presence of such bulk spinors is practically unavoidable, since they arise as 
superpartners of moduli fields. We assume the existence of an arbitrary number, N, of 
bulk Dirac spinors Q ^/ where I,J,--- = 1, • • • ,N. Given these fields, one might want 
to write the most general five-dimensional action bilinear in the spinors and respecting 
five-dimensional Lorentz-invariance. However, there is an important generalization to this 
that arises from including bulk gauge fields under which the fermions are charged. 
Within string- or M-theory such bulk vector fields arise quite generically as part of five- 
dimensional vector super-multiplets. Also gauging of the bulk fermions with respect to 
these vector fields is realized in large classes of models such as the five-dimensional brane- 
world models originating from heterotic M-theory [Q, ^ . Later on, we will compactify the 

■^Our conventions for five-dimensional spinors are summarized in Appendix 
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five-dimensional theory on a certain vacuum state to four dimensions. All we really have 
to require is that this vacuum state respects four-dimensional Lorentz invariance. As a 
consequence, in such a vacuum state, non-vanishing vacuum expectation values (VEVs) 
for the components ^4 of the bulk vector fields in the direction to be compactified are 
allowed. Hence, we should consider two types of bulk mass terms. First, there are mass 
terms that respect five-dimensional Lorentz invariance. We can think of such mass terms 
as being generated by VEVs of bulk scalars. Secondly, there are mass terms generated 
by VEVs of the components ^4 of bulk vector fields. They arise once five-dimensional 
Lorentz invariance is spontaneously broken to four-dimensional Lorentz invariance by the 
vacuum state. We split the five-dimensional action as 

— "Sbulk + "S'brane (2-1) 

into a bulk and a brane part. Taking into account the two types of mass terms as discussed, 
the relevant part of the bulk action is given by 



•Sbuik = d X 



^n'^ida^i + /xfj^-z^-j - ^{Mfj^j^j + h.c.) 



(2.2) 



Here, and Ms are the Lorentz-invariant Dirac and Majorana mass matrices, respec- 
tively. The subscripts "5" refer to their possible origin as being generated by VEVs of 
bulk scalar fields. These mass matrices have two vector-like counterparts, and My, 
generated by VEVs of bulk vector fields in the way indicated above 0. To see in some 
more detail how this happens, we note that in five dimensions the spinors ^ and their 
charge conjugate transform in the same way under the five-dimensional Lorentz group. 
Therefore, in addition to the "ordinary" kinetic terms of the form ^'y°'ida^, we can also 
write Majorana kinetic terms of the form ^^^y^ida^ + h.c. In the most general case one 
should gauge a linear combination of these terms. One can see that with an appropriate 
redefinition in the (^', space the Majorana kinetic term can always be removed. This 
is why, in the above action, we have only written the standard kinetic term. However, in 
general, one cannot simultaneously convert the couplings to the vector field to the same 
standard form. In other words, once the kinetic term is in the canonical form, the gauge 
transformations can still mix the ^ and ^''^ fields Q Therefore, after taking into account 
the vector field VEVs, we indeed generate the Dirac and the Majorana vector mass terms 



given in (2.2). From the structure of the various terms in ( p.2| ) we deduce the following 



symmetry properties of the associated mass matrices. 

^^S = lis^ , Ms = Ms^ , /iy = /iy^ , My = My^ (2.3) 

Are there other sensible generalizations of the bulk action that we have not yet taken 
into account? In models where the fifth dimension is divided by an orbifolding symme- 
try one should impose invariance of the action only under those infinitesimal coordinate 



'^In the following, we will use "S" and "V" as upper or lower indices depending on convenience. 

''As an example, one can consider the case of one bulk family and a U(l) gauge group acting on the 



fields by rotating the vector (>!', '5'^). 
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transformations that are compatible with the orbifolding. Since this may be a smaUer 
set of transformations it leads to generalizations of the action given above. Typically, 
in these cases, one can allow explicit orbifold-dependent functions such as step-functions 
in the bulk action. Another way of generating mass terms which effectively depend on 
the additional dimension is to consider a vacuum state with a non-flat metric. Explicit 
coordinate dependence of the mass terms generated in these ways might have interest- 
ing consequences for the structure of neutrino masses [28|. For example, generically one 



expects that the Kaluza-Klein part of the mass matrix does not have a block-diagonal 
structure any more. In this paper, we will not consider this interesting possibility further 
but we remark that the general formalism for the diagonalization of the mass matrix in 
Appendix ^ is equipped to handle even the more general case of a non block-diagonal 
mass matrix. 

Restricting to the case of coordinate-independent mass parameters, the bulk Lagran- 
gian ( p.2[ ) constitutes the most general expression compatible with four-dimensional Lorentz 
invariance and, therefore, contains all possible mass terms. However, we can ask whether 
we can generate these mass terms by mechanisms different from the ones already men- 
tioned. To do this, we note that the spinors ^/ transform in the fundamental of Sp{4:) ~ 
S0{5). A spinor bilinear, therefore, transforms under Sp{4) as 

4x4 = 1 + 5 + 10. (2.4) 

We have already discussed the coupling to the singlet 1 and the vector 5 leading to the 
scalar and vector masses, respectively. In addition, there can be a coupling to a second- 
rank antisymmetric tensor field transforming in the representation 10. However, only the 
components of this antisymmetric tensor field pointing into the internal space can take 
VEVs if we are to preserve four-dimensional Lorentz invariance. Clearly, for the case of 
one additional dimension these components must vanish by antisymmetry of the tensor 
field. Therefore, mass terms are generated by VEVs of scalar and vector fields only. 

Next, we should set up the brane action. We choose x^, where = 0,1,2,3 

as the coordinates longitudinal to the brane and ?/ = as the transverse coordinate. 
Furthermore, we locate the brane at y = 0. The essential fields on the brane that we need 
to consider for our purpose are the left-handed lepton doublets Lj, where z = e,^,r and 
the relevant Higgs field H. As mentioned in the introduction, we will not consider singlet 
neutrino fields on the brane. We write the lepton doublets in a basis in family space where 
the charged lepton mass matrix is diagonal and positive. The most general brane-bulk 
Yukawa couplings invariant under the standard model group can then be written as 



'brane 

'{y=0} 



/ 

J>.v 



d^x 



/M5 



"^iLiH - -ji^'^jLiH + h.c. 



(2.5) 



Since the bulk spinors with the normalization fixed by the bulk action ( |2.2| ) have mass 
dimension two, these couplings are suppressed by the square root of a mass scale M5 
defined by 

2,RM, = (2.6) 
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where M^, is the string scale. Furthermore, R is the radius of the circle on which we 
compactify the fifth dimension transverse to the brane in preparation for our reduction 
to four dimensions. Also, we restrict the range of the corresponding coordinate as y G 
[—ttR, ttR] with the endpoints identified. 

As mentioned in the introduction, in the context of string- or M-theory, brane-bulk 
couplings may arise from higher-dimensional gravity-induced operators after inserting ap- 
propriate VEVs. From a theoretical viewpoint, their magnitude is rather uncertain and 
is, in our case, parameterized by the dimensionless couplings hji and h'f-. In the context 
of string theory, we have, of course in mind that our five-dimensional model is obtained by 
compactification of the ten (or eleven) dimensional theory on a five- (or six-) dimensional 
internal space with radii of compactification smaller than R. In the simple case in which 
all those internal radii are equal and given by p, the quantities R, p and the string scale 
M* are related by 

167r(27ri2M,)(27rpM,)5 = (^^^ . (2.7) 

From this relation, in models with a low string scale M*, we will need p~^ to be smaller 
than the string scale in order to get the correct value for the low-energy Planck scale. In 
such cases, the five-dimensional action specified by eq. (|2.1| ) should be considered as an 
effective theory below p~^ (rather than below M*). This also means ^ that we should 
restrict bulk masses to be smaller than p~^. 

If the bulk fermions propagate in all ten dimensions, as is suggested by their possible 
interpretation as modulini, then, from eq. (p.6|), the dimensionless couplings hji and hj- 
are simply given by their ten-dimensional counterparts. On the other hand, if the bulk 
fermions propagate only in 5 of the ten space-time dimensions (where 5 < 5 < 10) the 
couplings h are given by /i = /i5(27r/>M,)(i°-^)/2. 



2.2 The four-dimensional effective action 



We would now like to reduce the above five-dimensional theory and calculate the 
associated effective four-dimensional action. We require the bulk fields to be periodic on 
the circle ^, that is 

^i{-itR) = ^i{itR) . (2.8) 
Correspondingly, the bulk fermions can be expanded into Kaluza-Klein modes as 



^We would like to thank R. Barbieri for reminding us about this. 

® Other possible choices are obtained by imposing generalized periodicity conditions of Scherk-Schwarz 
type Those will lead to complications in the structure of four-dimensional masses for the bulk fields. 
Particularly, there can be a shift in the Kaluza-Klein masses. However, we will not consider these exotic 
possibilities here. 
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In addition, it is useful to express each of the modes ^ni in terms of two four-dimensional 
(left-handed) Weyl spinors C-ni, Vni as 



nl 



Vnl 



(2.10) 



Here the conjugate ^ of a Weyl spinor ^ is defined as ^ = e^* where e is the two-dimensional 
epsilon symbol. For further details of our conventions we refer to Appendix^. The four- 
dimensional effective action associated to the brane- world theory (p.l|), (2.2), (2.5) is then 
given by 



where the kinetic terms of the bulk fields read 

^kin = ^ [vii^^^id^Vni + iijcr^id^ini 



(2.11) 



(2.12) 



The all-important mass part of the Lagrangian can be conveniently written in the form 

1 



E 



Wn Mr 



rj-n 



+ h.c. 



(2.13) 



where, for ease of notation, we have adopted a vector notation in family space. That is, 
for example, rjn = (??ni, • • • , VriN)'^ and u = {ue, i'^, i't)'^- Given our conventions for Weyl 
mass terms in Appendix |A|, here, the transpose T refers to family space only. Written in 
block-form, the mass matrix 7W„ for the Kaluza-Klein modes at level n takes the form 



Mr 



R 



Ma 



(2.14) 



Here Mi and M2 are Majorana mass-matrices related to the scalar and vector Majorana 
mass-matrices appearing in the five-dimensional bulk action (|2.2D by 



Ml 



-Ms + My , M2 = {Ms + MvY 



(2.15) 



The Dirac mass-matrix /j, encodes both scalar and vector Dirac masses and is explicitly 
defined by 



(2.16) 



From eq. ( |2.3| ), it is clear that Mi and M2 are symmetric matrices while, generally, there 
is no further constraint on /x. 

Note that the bulk mass matrices satisfy the relations 



M. 



Mi 



(2.17) 



They ensure that the Majorana mass matrix 



M 







'10 so 


'l—n s— n 


'In Sn 


V 
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••• 






Vo 
Co 


N 


Mo ■■■ 










V-n 
C— n 


N 








■' '■n 


Vn 


N 





Mn 






(2.18) 



associated with the mass Lagrangian (2.11) is indeed symmetric as it should be. 
The brane-bulk couphng matrix N has the form 



m 
m 



(2.19) 



The mass matrices m and mP are related to the brane-bulk couplings in the brane ac- 
tion (|2.5|), respectively and are explicitly given by 



mii 



vhii 



V27ri?M5 



hiiv 



Mpi 



li 



li 



V2nRM5 



h'jiV 



Mpi 



(2.20) 



Here is the Higgs vacuum expectation value. 

In what follows we will have occasion to choose the brane-bulk couplings, hji, to be 
of a given size for phenomenological reasons. However, there are lower bounds on M^/hu 
at tree level or on Mjh\ at one-loop-level by requiring that there be no significant 
departure from Standard Model processes [^]. It turns out that these constraints can be 
satisfied within the explicit model to be discussed later. Moreover, it has been suggested 
that values of hn larger than one may be inconsistent with staying in the perturbative 
domainj^. However, from eq. (2.7), we see that the string scale is very sensitive to the 
nature of compactification from higher dimensions so it is always possible to arrange for 
hii to remain perturbative through a choice of M^, . 



2.3 Symmetries of the brane-world action 

For the purpose of model building it is important to consider the various symmetries 
that can be imposed on the general brane-world action ( |2.1| ) in order to forbid some of the 
mass terms. A general classification of these symmetries can be obtained by analyzing the 
invariance properties of the bulk kinetic terms. Instead of going through such a systematic 
procedure, here we merely list some of the most important symmetries. The most obvious 
candidate is the part of the five-dimensional Lorentz group connected to the group identity 

'''We thank J. Valle, A. loannisian and A. Pilaftsis for comments on this point. 
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which we call Cq. Infinitesimal transformations of this type act on the fields in our model 

as 



5^1 = '-uj'^^^^p^ii , 8U = ^LO^'^f.uLi (2.21) 

where uj°'^ are infinitesimal parameters and are the generators defined in Appendix^. 
The original five-dimensional action should, of course, be invariant under this symmetry. 
However, as discussed above, this is not necessarily the case for the vacuum state of the 
five-dimensional theory. Imposing this symmetry, therefore, means requiring a vacuum 
state which is invariant under infinitesimal five-dimensional Lorentz transformations. Cor- 
respondingly, this symmetry forbids the two vector-like bulk mass terms in the action ( |2.2D 
that we have attributed to spontaneous Lorentz-symmetry breaking by the vacuum state. 
On the other hand, it allows all other terms. 

Furthermore, we can consider the generators of those parts of the five-dimensional 
Lorentz group that are not connected to the identity. An important example is the parity 
transformation 

P5 : y^-y (2.22) 

in the direction transverse to the brane. On the fields it acts as 

^/(y) ^ 75^K-y) , U^Li. (2.23) 

In principle, the transformation law for has a sign ambiguity. However, and its 
charge conjugate ^'^ transform with opposite signs. We can therefore remove this sign 
ambiguity by exchanging and appropriately and calling those fields that trans- 
form with a positive sign under P5 as in eq. ( p. 23 ). It is also useful to have the action of 



this symmetry on the Weyl Kaluza-Klein states appearing in the four-dimensional effective 
action available. It is given by 

Cnl i-nl , Vnl ^ -rj-nl , Vi ^ Ui . (2.24) 

This Z2 symmetry forbids bulk Dirac mass terms and the brane-bulk couplings rrf. It 
allows all other mass terms. 

For N bulk spinors the bulk kinetic terms have a global U{N) symmetry. Particularly 
interesting for model building are the U{\) subgroups acting on the spinors as 

$7 ^ e'^'^'^i , Li e*"'''L, (2.25) 

where a is a continuous parameter and Qi and qi are charges. The group action on the 
Kaluza-Klein modes reads 

Cn/ ^ e-*"«^en/ , ry„/^e'"%n/, Vi^e^^'i'vi. (2.26) 

A specific example is lepton number defined hy Qi = qi = 1. This symmetry allows 
the Dirac mass terms and forbids the Majorana mass terms. Clearly there are many more 
interesting choices that can be made in order to constrain the mass terms in a sensible way. 
We will encounter another example later on. In the context of string- or M-theory, if one 
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-'-0 






1 I A • ,S' 

scalar Dirac ^jj 


inv. 


not inv. 


Qi = Qj 


scalar Majorana M^j 


inv. 


inv. 


Qi = -Qj 


vector Dirac /i^j 


not inv. 


not inv. 


Qi = Qj 


vector Majorana Mjj 


not inv. 


inv. 


Qi = -Qj 


brane-bulk m/j 


inv. 


inv. 


Qi = Qi 


brane-bulk mfj- 


inv. 


not inv. 


Qi = -Qi 



Table 1: Invariance properties of the mass terms with respect to the symmetries introduced 
in the text. For the U{1) symmetries the conditions for invariance are stated. 



thinks about these symmetries as being global, it appears more likely to have invariance 
under a discrete subgroup of U{1) rather than under the fuh continuous symmetry. Hence, 
when we consider U{1) symmetries in the following, we will always have in mind that one 
might have to replace those symmetries by appropriate discrete subgroups. 

In Table || we have listed the various mass terms and their transformation properties 
under the above symmetries in order to give an overview of the constraints available for 
model-building. 

A final remark concerns models that are obtained by compactification on the orbifold 
ITL2 where the action of the Z2 symmetry is specified by P5. While, so far, we have 
just considered the reduction on a circle, the orbifold case does, in fact, not need to be 
treated separately for our purposes. This is for the following reason. We can consider 
our general four-dimensional model specified by the Lagrangian ( p. 11 ) and impose P5 as 
a symmetry. Using the P5 transformation properties ( |2.24| ) of the four-dimensional fields, 
we can then define Z2 even and odd modes. Our Lagrangian ( |2.11 ), being bilinear in the 
fields, then decomposes into two parts, one which involves only even fields the other one 
involving only odd fields. The ordinary neutrinos Vi are even fields and, therefore, only 
appear in the even part of the Lagrangian. As a consequence, the odd modes decouple 
from the physics that we are considering in this paper. Orbifolding such a Z2 invariant 
model merely means that the already decoupled odd states are truncated. We therefore see 
that for each orbifold model there is a corresponding P5 invariant model with equivalent 
physics as far as neutrino masses and mixing are concerned. 



3 General structure of masses and mixing 

It is clearly desirable to obtain general results about the masses and mixing in our 
brane-world model that are independent of specific model-building considerations. We 
would therefore like to present a number of general conclusions derived from the effective 
four-dimensional mass Lagrangian ( |2.13| ) before discussing explicit models. 



3.1 Masses of bulk modes 



As a first step, we will analyze the bulk part of the mass Lagrangian separately. This 
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bulk part is given by the first term in eq. ( 2.13 ) and reads 

-^bulk = 



^n)Mn(^^2^+h.C. (3.1) 



where the mass matrices have been defined in eq. ( 2.14| ). Clearly, the mass diag- 



onalization of this bulk action will have to be corrected when the brane-bulk couplings 
are taken into account. However, usually we will be interested in the case where these 



couplings and, hence, the corresponding corrections are small. In this case, studying (3.1), 
provides us with useful information about the approximate spectrum of bulk states. 

Let us first note that the mass matrices can be diagonalized by unitary matrices 
where n G Z, such that 

Mn = U^MnU.n (3.2) 

are diagonal matrices with entries 

TWn = diag(M„i+,--- ,MnN+,Mni-,--- ,MnN-) ■ (3.3) 
The original states {r]n,(,n) are related to the mass eigenstates {fjm^n) by 

The eigenvalues and the matrices C/_„ (up to unitary rotations within each eigenspace) 
can be computed by diagonalizing the hermitian matrices AinA4n- 

For each n > we have AN Weyl mass eigenstates {f]n-,'f]-n-,in-,i-n)- It is immediately 



clear from the structure of (3.1) that the components with opposite mode number n 
arrange themselves into Dirac states. For each integer n > 0, we, therefore, have 2N 
Dirac states with masses denoted by {Mni+-,Mni-)i=i^--- ,n- This result was, perhaps, not 
quite so obvious from the original five-dimensional bulk action ( p. 2] ). We remark that the 
corrections induced by the brane-bulk couplings (to be included below) will generally not 
respect this Dirac structure. However, for small couplings we will still retain approximate 
Dirac states. The modes for n = 0, on the other hand, do not generally group themselves 
into Dirac states and represent 2A'" Weyl states with Majorana-Weyl masses denoted by 
(Mo/+,Mo/_)/=i,.., TV- 
It is difficult to extract more information about the spectrum of bulk states without 
further assumptions. To see the range of possibilities, we will, therefore, consider two 
specialized cases. First, we consider a situation where the vector mass matrices are zero, 
that is jiy = My = 0, and the scalar mass matrices ^5, Ms are real. Then we have 



m 



Hence, the matrix Mq for the lowest mode n = is symmetric and can be diagonalized 
by an orthogonal matrix V leading to A^o = diag(Mo/+, Mo/_)/=i,... = V'^M-qV. Then 
one can easily see that the matrices Un defined by 

V for n = 
Un= { PV for n > (3.6) 
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with 



J_( Itv liv 

^ V liV 



(3.7) 



diagonaUze the mass matrices 7W„ for arbitrary mode number n. From eqs. (|3.5| ) we then 
find for the eigenvalues (neglecting phase factors) 



Mni± = \jMl^+'^. (3.8) 

We recall that, in this equation, Mqi± are the masses of the n = modes which depend 
on the scalar Dirac and Majorana mass matrices. We can therefore express the complete 
mass spectrum as well as the diagonalizing unitary matrices in terms of the corresponding 
quantities for the n = mode. From the general discussion above it is clear that the 2N 
masses given by eq. (|3.8| ) correspond to 2N Weyl states for n = and to 2N Dirac states 
for n > 0. The spectrum (|3.8[) has the form one usually expects for a Kaluza-Klein tower 
of particles. In particular, the masses Mqi± and, hence, the scalar Dirac and Majorana 
masses, provide a lower cutoff for the spectrum. Therefore, in order to have light bulk 
masses, the scalar mass terms must be small compared to the string scale. As discussed 
earlier, this is the natural expectation if the bulk fermions are interpreted as superpartners 
of moduli fields. As we will see below, the situation is quite different for the vector-like 
masses as they do not represent lower bounds on the bulk mass spectrum. 

In our second example, we would like to include all types of mass terms. In particular, 
we are interested in the effect of the vector-like masses which we have previously set to 
zero. As a simplification, we restrict the analysis to the case of one family, that is = 1 
and assume that all mass parameters are real. We are, therefore, dealing with 2x2 mass 
matrices 



f -Ms + Mv fis- W - 1 A 



(3.9) 



In this case, the two eigenvalues of AiliM-n are given by 



M2^ = M| + /i| + M2 +^2^ + ^±2^/i2.(^^ + M|) +M2 (^^ + M| + /i| 

(3.10) 

The first five terms on the left-hand side of this equation are as expected and are of the 
form ( |3.8| ) that we have encountered in the absence of vector-like masses. The square root 
term (which vanishes for vanishing vector-like masses) does not fit into this structure. 
Specifically the eigenvalues corresponding to the negative sign in front of the square root 
are not bounded from below by the vector-like masses. To see this in more detail let us 
consider the four cases in which only one type of mass term is present. For a non- vanishing 
scalar Dirac mass but all other mass terms vanishing we have 

^n±=4+j- (3.11) 
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The spectrum is of the form eq. (|3.8| ) and is indeed bounded from below by ^s- Similarly, 
this is the case for a scalar Major ana mass and all other mass terms vanishing. In this 
case we have 

Ml^ = M| + ^ . (3.12) 

The situation is quite different, however, for the vector-like mass terms. Let us consider a 
situation with a vector Dirac mass only. Then the spectrum is given by 



M4= (^|^y|±i^J . (3.13) 
Analogously, for a vector Majorana mass as the only mass term we have 



Ml^=[\Mv\±^-'^ . (3.14) 

Hence, in these last two cases, the spectrum is not bounded from below by the vector-like 
bulk mass term. While the mass of the n = mode is still given by the vector bulk mass 
(as in the scalar case) the masses of some of the higher modes are lower. In fact, the 
modes n specified by the integer part of ±fisR i^MsR) have the smallest mass given by 
US modulo {Ms modulo R~^). As a consequence, there will always be modes with 
masses of the order 1/R independent on the size of the vector bulk masses. 

3.2 Including brane-bulk couplings 

So far, we have analyzed the spectrum of the bulk modes for vanishing brane-bulk 



coupling. However, our main task is to diagonalize the full mass matrix (2. IS). A general 
formalism for an exact diagonalization is set up in Appendix |^. However, in many cases 
it is sufficient to use a perturbative diagonalization leading to a see-saw type mechanism. 
This allows one to bring the matrix ( |2.18| ) into block-diagonal form. The approximation 
is also well-suited to understand the physical implications of our models intuitively. To 
see how this works, it is useful to introduce the quantities 

En = M-^M . (3.15) 

For a see-saw diagonalization to be sensible these quantities should be small compared to 
one. However, we should also take into account that we have to sum over a large number 
of modes n. Therefore, the more precise requirement is that 

e' = Y.eien (3.16) 
nez 

should be small compared to one. Provided this is the case, we can perturbatively diago- 
nalize the mass matrix neglecting terms of order and higher. To do this, we define the 
"see-saw" basis {v' ,£,'^,rj'^) by 

(3.16a) 

nGZ ^^"^ 
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Diagonalizing, the mass Lagrangian (|2.13| ) then takes the block-diagonal form 



>Cmass = -\u'^m'u' - \ [riJ Mn (fA . (3.17) 



Here, the mass matrix ml for the light neutrinos is given by 



I 

m 



Y,^M-'N. (3.18) 



These results hold to order e and receive corrections at order e^. In a second step we 
now need to diagonalize the various blocks in ( 3.17] ). To do this, we introduce the mass 



eigenstate basis {v,'r)n:^n) related to the see-saw basis by unitary rotations 



such that the mass matrices 

rh = U^m'U , Mn = U^MnU^n (3.20) 

are diagonal. The diagonalization of the bulk states is exactly as discussed in the previous 
Subsection and the unitary matrices [/„ are precisely as introduced in eq. (|3.2D . We can, 
therefore, use the results for the bulk mass spectrum and the matrices Un obtained there. 
The diagonalization of the light mass matrix m' is clearly more complicated and will, for 
particular cases, be explicitly carried out later on. With these quantities available we can 
now express the electroweak eigenstates u in terms of the mass eigenstates {v-,fjniin) as 



U{^ + Y.'^-nUn(f] (3.21) 



where we have neglected terms of order e^. This equation can, of course, also be derived by 
expanding the exact formula (B.13) to first order in e. The result (3.21) is, in fact, central 



to the subsequent discussion of the oscillation phenomenology. Provided that the quantity 
defined in eq. (3.16) is much smaller than one, as we have assumed, this equation asserts 



that the electroweak eigenstates are predominantly composed of light mass eigenstates 
z) but have small admixtures of the heavy mass eigenstates {fin,£,n)- The size of these 
admixtures is determined by the quantities e„. The mixing matrix U then plays the role 
of the standard (MNS) neutrino mixing matrix. This offers the interesting possibility that 
some of the observed oscillation phenomena could be attributed to oscillations into the 
light states (playing the role of the standard neutrinos) while other phenomena could be 
due to oscillations into the heavy states (which are predominantly bulk states). In fact, 
given the various mass terms in our model, the parameters in these two sectors can be 
largely independent. 

To verify this last statement, it is useful to evaluate some of the above quantities more 
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explicitly. We introduce the matrices 
Xn = 



'r 



Ml U + 



'r 



Mo 



-1 



m 

'r 



Rn — 
Sn = XnMi (^/i + 

and further note that R-n = R^ and S-n 
mixing parameters e„ of the bulk states 



(3.22) 
(3.23) 

1 

(3.24) 

ST . With these definitions we find for the 



M2X„ 



-RnVn'^ + Xt^m 
X„m^ — Srifn 



In addition, the mass matrix m' for the light states is given by 



m 



' = ^ \rrf^ Rnirf + Snm — m!^^ m — Xn'm'^ 



(3.25) 



(3.26) 



Unfortunately, it is difficult to explicitly sum this series in the most general case. However, 
partial results for special cases will be given below. We would now like to specialize these 
results to models that preserve standard lepton number |20, 22] with charges Qj = qi = 1. 
From Table |, we see that in this case all Majorana-type mass terms, that is Ms, My and 
m'^ must vanish. As a consequence, we have Rn = Sn = and m' = 0. Therefore, we have 
three massless Majorana states. Being a direct consequence of imposing the lepton number 
symmetry these states are exactly massless even beyond the see-saw approximation that 
we are considering here. The mixing parameters take the simple form 







m 



(3.27) 



Prom eq. ( |3.21 ) we, therefore, see that the electroweak eigenstates are generally non-trivial 
superpositions of the massless states and the bulk states. The mixing matrix U in eq. ( p. 21 ) 
is, of course, arbitrary and corresponds to a redefinition of the degenerate massless states. 
The size of the bulk components n in this superposition is controlled by the quantity 
(/i + in/ R)~^m. For Dirac masses jj, ^ m these bulk components are small and the 
electroweak eigenstates are predominantly given by the massless states. This has to be 
contrasted with the opposite case /i ^ m or rather = that has been mostly considered 
in the literature so far. In this case, the massless states are simply the zero modes of the 
Kaluza-Klein towers and they decouple from the electroweak eigenstates. Moreover, the 
masses of the lightest modes are then given by the brane-bulk mixing m. In contrast, in 
the presence of a bulk Dirac term, the masses of the lightest bulk modes are determined 
by this Dirac mass. Keeping the Dirac mass fixed, we can then vary the brane-bulk 
coupling m thereby changing the importance of the heavy components in the electroweak 
eigenstates and, at the same time, keeping three massless, non-decoupled states. Hence, 
the lepton-number conserving mass terms control the oscillation phenomenology associated 
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with the bulk states. Small lepton-number breaking Majorana-type terms, on the other 
hand, generate small masses for the formerly massless states and fix the mixing matrix 
U in eq. ( 3.2l| ). Therefore, they control the oscillation phenomena related to the light 
states. This confirms the picture in which there are two different sectors, both potentially 
relevant for oscillation phenomenology, and which, by choosing the various parameters in 
the model, can be controlled quite independently. As we will see, to some extend this 
is also the case for the model with bulk Majorana terms to be considered in Section |^ 
although there the Majorana term determines the mass of the lightest bulk states and is 
also involved in the expression for the light neutrinos. 



3.3 Model building options and scales 

We have seen that our models typically accommodate two sectors of mass eigenstates, 
both of which generically couple to the electroweak eigenstates. In particular, there are 
three light eigenstates which (in the see-saw approximation) are comparable to the mass 
eigenstates in standard three-neutrino models. In addition, we have Kaluza-Klein towers 
of heavy states. Oscillations into both sectors may play a role in explaining the observed 
oscillation phenomena. A broad classification of the various possibilities can be given using 
the relevance of the Kaluza-Klein states for oscillation phenomenology as a criterion. The 
relevance of the Kaluza-Klein sector is basically controlled by the compactification scale 
1 /R and the smallest mass Mmin in the bulk Kaluza-Klein spectrum |^. These should be 
compared to the mass differences. Am, needed to explain the various oscillation phenom- 
ena. Normally, then, we take Am to be in the range relevant for either the MSW solution 
of the solar neutrino problem, that is, Am = 10~^-10~^eV or in the range for oscillations 
of atmospheric neutrinos, that is. Am = 10~^-10~^eV. We can now distinguish three 
cases. 



1/R ^ Am and M^[^ > Am : In this case, none of the Kaluza-Klein states play 
a direct role for the observed oscillation phenomena. Therefore, these phenomena 
have to be explained entirely by oscillations between the standard neutrinos, that is, 
by the first term in eq. (3.21). In such a situation, it may seem that the brane-world 
nature of the model does not play any role at all. However, as can be seen from 
eq. ( p. 26 ), these states do affect the mass matrix of the light neutrinos and, hence, 
the standard mixing matrix U. As we will discuss in more detail below, we can 
have models which mimic ordinary see-saw models with the mass scale of the right- 
handed neutrinos set by a bulk mass term or models with "geometrical" see-saw 
mechanism [10| where the right-handed mass scale is given by 1/R. In summary, 
this case may, therefore, lead to interesting alternatives to conventional "see-saw" 
models of neutrino masses. 



• 1/R ^ Am and Mmin < Am : In this case, only the bulk modes with the mini- 
mal mass A/jYiin are directly accessible for oscillations. Such models are, therefore, 
equivalent to conventional models containing a (small) number of sterile neutrinos. 

*In cases with vanishing vector-hke bulk masses this smallest mass coincides with the mass of the n — 
mode. For non-v anishing vector-like masses, however, this is not necessarily the case. See the discussion 
around eq. (|3.1C|) for further details. 
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We remark that an explanation of the atmospheric neutrino deficit due to oscilla- 
tion into a sterile neutrino appears to be disfavored ||2^. Given this one must try 
to use this case to explain the solar neutrino deficit by sterile neutrinos and the 
atmospheric deficit (and, perhaps, the LSND result) by standard oscillations. How- 
ever, large mixing angles between active and sterile neutrinos might be problematic 
from the viewpoint of supernova physics as well as standard model universality con- 
straints ||l^. On the other hand, the small mixing angle solution for oscillation into 



a sterile neutrino seems to be disfavoured by recent Super-Kamiokande data |30]. 
Hence, this case, although perhaps not excluded does not appear to be very attrac- 
tive. 



• 1/R< Am and Mmin ^ Am : In this number of bulk modes can contribute 

to oscillation phenomena. The case that Ve or has a large bulk component in order 
to realize the large angle MSW solution is problematic for the same reasons as above. 
On the other hand if the mixing of bulk modes with the electroweak eigenstates is 
small then the small-angle MSW solution of the solar neutrino problem might be 
attributed to oscillations into bulk states. As we will discuss, the phenomenology 
is somewhat different from the second case and, in fact, compatible with the recent 
Super-Kamiokande results as well as with the other solar neutrino experiments. 
At the same time, the other oscillation phenomena should then be explained by 
oscillations between the standard neutrinos. 

Some of these options will be discussed much more explicitly in the context of the 
specific model to be presented in the following Section. 

What about the choice of the string scale M5, the radius R of the fifth dimension and 
the radius p of the other five (or six) dimensions that have already been compactified to 
obtain the five-dimensional model? In models without any explicit bulk masses these three 
parameters are basically fixed if one of the oscillation phenomena is to be explained by bulk 
states because one has to obtain the correct value of the low-energy Planck mass as well 
as appropriate masses and mixing of bulk states. In the presence of bulk masses, however, 
things are less constrained. For example, as we have seen from the above discussion, the 
value of 1/R is not fixed and merely determines how many bulk states actively contribute 
to the oscillation. Thus, the phenomenological possibilities are much richer when bulk 
masses are allowed. 



4 Supernova 1987 a 



One may obtain significant bounds on the mixing angle of a sterile neutrino from the 
condition that sterile neutrino emission should not cool the supernova too much, that is, 
sterile neutrinos should not carry away more energy than is carried by doublet neutrinos, 
approximately 10^^ ergs. Production of sterile neutrinos proceeds by either incoherent 
processes or through coherent oscillation of the active neutrino state into the sterile Kaluza 
Klein tower. We shall discuss both these processes. As we shall see the latter places the 
most stringent bounds on the properties of the sterile neutrino. 

As a preliminary let us summarise the salient points [31|. The core of the super- 
nova of approximately 10 km radius reaches nuclear densities with central density pQ ~ 
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3.10^^ gcm~^. The shockwave moving into the core produces a hot gas of neutrons, pro- 
tons, leptons and radiation. Due to the high density neutrinos are trapped within the core 
and following trapping and the establishment of chemical equilibrium, the lepton fraction, 
Yl, is given by Y;, ~ 0.35. The relative contribution to of e and Ve is determined 
by the equilibrium condition on the chemical potential. Initially the lepton fractions are 
constant throughout the core with Yg = 0.28 and Y^^ = 0.07, corresponding to Fermi 
energies /Ug ~ 235 MeV and fi^e ~ 180 MeV. Thermal distributions of electron neutrinos 
are maintained by the dominant charge exchange process 

ep ^ V(.n . (4.1) 

The remaining species of doublet neutrinos are produced by the processes p2| , ^] 

e+e" — > vv (4.2) 
nn — > nnvv . 

We will assume that the mixing angle of the sterile neutrino with the active neutrino 
is small enough that the sterile neutrino can escape from the core. The case of larger 
mixing angles has been studied in References |34, 35| and leads to a lower bound > 0.02 



for the square of the mixing angle in the medium to the part of the KK tower with the 
appropriate mass difference. Apart from an isolated state which may be near resonance, 
this corresponds to much larger mixing angle in vacua, inconsistent with universality 



constraints at least for the electron and muon neutrinos |15|. 

In calculating the sterile neutrino production rates we need to distinguish between v^ 
and Vt production for which the chemical potential is small ^ and z/g which has a large 
chemical potential. 

4.1 Incoherent production from z/„ and v^-- 



We first consider the incoherent production of sterile neutrinos from fi and r neutrinos |34, 
p!9| . Incoherent production takes place if the SM neutrinos are Dirac, with the right- 
handed component being sterile. We note that in the model discussed in the next Section, 
the SM neutrinos are on the contrary mainly Major ana- type and so the bound is much 
weaker. Since the energy stored in the thermal distribution of fi and r neutrinos is a 
very small proportion of the overall energy in the core it is important to first determine 
the processes capable of transferring energy from the electron and nucleon sector to the 



neutrino sector. The dominant production processes are those in eq. (O). It turns out 
that these processes give similar rates so here we will content ourselves with a simplified 
discussion assuming production via the electron positron channel. The electrons and 
positrons are in thermal equilibrium with a gamma radiation which in turn is in equilibrium 
with the hot matter behind the shock, that is, with the neutrons. These processes are 
much faster than weak and thus we can consider them as instantaneous, maintaining a 
thermal distribution of electrons and positrons. Following reference the v rate of 
energy loss is given approximately by 

Qe = ^^((1 - Cy f + Cl)m6 (^)' ( nee-^=A^ (4.3) 



'•'As discussed in Ref. [p8| this may be too strong an assumption. 
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where G is the Fermi constant, m is the electron mass, re/ = Nf — Nj where Nf is the 
number density of fermions and (1 — Cy)^ + « 1/4. Clearly Qe is strongly peaked at 
high temperatures so neutrino production is dominated by the high initial temperature 
produced by the shock wave, T ~ 30 MeV. To obtain a useful bound it is necessary to 
discuss what happens to the neutrinos produced in this process. The production cross 
section Qe involves electrons with energy of 0{fie) while the positron energy integral is 
approximately proportional toj^"" dE E^ e~^/^. One may see that it predominantly comes 
from the energy range E > 3T. Thus the produced neutrinos have energy much greater 
than the average temperature. There are two processes relevant to the fate of these 
energetic neutrinos 



Elastic scattering from the neutrons. This process causes the neutrinos to change 
direction but does not significantly change their energy. It is the dominant process 
giving rise to neutrino trapping but not thermalisation. It also gives the most efficient 
process for incoherent production of sterile neutrinos. Following references 34] 



this happens at the rate {m/E) Tnnc where Tnnc is the neutral, current collision 
rate on free nucleons 

Tnnc ^ -GINbEI . (4.4) 

TT 

and m is the brane bulk mass, as defined in eq. (2.20). 



Elastic neutral current scattering processes of the neutrinos off electrons, positrons 
and neutrinos or neutrino annihilation processes. The most rapid process is the 
scattering off electrons because of their greater abundance due to the large chemical 
potential. It will rapidly thermalise the neutrino, reducing the dominant elastic 



scattering of eq. (4.4). However the neutrino remains until it annihilates and during 
this time the production of sterile neutrinos continues. The rates for annihilation is 
given by 

— 

Tann ^ ^N,{T)EI (4.5) 

vr 



so 



~ 3.5 • 10^ ( j . (4.6) 



r Tnnc Nb q /30 MeV\ 

i -L a.nrt. ^ * I 



V 



Then, for small {m/E), the probability the neutrinos oscillate into sterile neutrinos is 
simply given byP = = (^)^ RE where the factor RE corresponds to the number of Kaluza 
Klein states energetically available. If we require that the energy loss QeVtP should be 
less than the energy carried off by doublet neutrinos we immediately obtain the bound 

< 3 . 1053 erg LEM[. ~ 5 . lQ3^^^eil/kT-l/k so MeV) (4 7) 

eV F QeVt 

where we have assumed a constant density of neutrons and chemical potential of electron 
neutrinos {fie — 200 MeV) over the volume V cc r^^^^ ~ (lO^cm)^ . Since the temperature 
in the supernova falls rapidly over the first ten seconds and the production cross section is 
very sensitive to the temperature, the dominant production occurs within the first second, 
that is for t = 1. In eq. (|4.7| ), T refers to the initial temperature after the shock and 
assumes it lasts for one second and we have taken E = 3T. 
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4.2 Coherent production of and i/r- 

A much more stringent bound follows from coherent oscillations of the active neutrino 
state into the Kaluza Klein tower of sterile neutrino states. There are two contributions 
that must be considered, transitions induced through level crossing and mixing-induced 



transitions. The former have been discussed in references [|19|, 12, 22]. Here we extend 
this discussion taking account of the non- vanishing chemical potentials for the neutrinos 
and their radial dependence. As we shall discuss there is an important limiting effect 
that reduces the effect of level-crossing conversion. The latter mixing-induced transitions 
are also much more important than the incoherent processes considered above. To see 
this note that the neutrino flavour state is a coherent mixture of active and sterile states 
with an average probability of finding the original flavour state Va in the sterile state 
given by sin^(20)/2 where 6 is the f^, Vs mixing angle. The mixing-induced transition 
rate is given by r„„c sin^(20)/2 to be compared with the incoherent production rate of 
^ nnci'fn / E)"^ . The mixing angle to the sterile neutrino tower in vacua is approximately 
m/Mn, much larger than the equivalent factor [m/E) relevant to incoherent production. 
In the medium of the supernova the mixing can be much larger due to resonant effects 
because the electron doublet neutrino can be degenerate with one of the KK levels. Thus 
the mixing-induced transitions will also be much larger than the non-resonant transitions 
and can compete with the level crossing contributions. Here we will compute only the 
level-crossing contribution; a detailed calculation of the full contribution will be presented 
elsewhere. 

Let us consider the muon neutrino case - an identical discussion applies to the tau 
neutrino. The muon neutrino effective potential, V^, in the supernova core is given by 

= V2G{2n^^ + + n^, - 0.5n„) . (4.8) 

Initially only n^^ and n„ are significant yielding a negative V^, of order —10 eV. In 
this case it is the antineutrino, 17^, that crosses the Kaluza-Klein levels through matter 
effects. It acquires a resonant mass m„ ~ n/R ~ {2Vf^Ey/^ ~ 10^'^ eV. The effect of 
these oscillations can be reliably estimated for small mixing angles, because the width 
of each resonance is smaller than the separation between resonances. Thus the survival 
probability for a standard neutrino produced in the core is given by the product of the 
survival probabilities in crossing each resonance P^u — ^nPn, where Pn — e~'^"'"^'^, and 

4m2 V ^ 
^ -Wlvjdr ■ ^'-'^ 

Following reference |jl^ we approximate the density profile by p{r) oc e~'~/^'=°'^<= . This gives 
V/{dV/dr) = Tcore and hence 7„ ~ m^/10~^ eV^. 

To determine the energy loss due to coherent processes we must determine the survival 
time of the active neutrino after production. The neutrinos initially produced by the 
annihilation processes eq. (^]^) thermalise through the elastic scattering e ^ f ^ e. The 
electron abundance is very large due to its large chemical potential and so thermalisation 
occurs on a short timescale compared to that for annihilation which stops the level- 
crossing production process. Thus most of the production of sterile neutrinos occurs after 
the (anti) neutrino acquires its thermal distribution. As the antineutrino moves towards 
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the surface of the supernova its mass decreases due to the density change and thus it 
crosses several of the Kaluza-Klein levels before it annihilates. To estimate the number of 
levels crossed we first compute the mean free path, X'{T,E^), before z^z7 annihilation. The 
annihilation cross section is 

a_ ^ _^^n,{E,,T)El (4.10) 



giving 



/>oo 

\'{T,E^)-^ = C / S|cT^^_^e+e- [ETj + E^)e-^^/'^dETj . (4.11) 
Jo 

It is convenient to eliminate the constant C by writing this as 

A'(r,£„)-' = KiT.E. = 100 MeV)- ■ { ,,^r,Z=mMeV)-^ ) ''"'^ 



( f°° E^(E- + E )'^e'^/^dE- 
A'o(r, E, = 100 MeV)-i - ' + 



E^e-^/^dE^ {100 UeVfJ 
where Aq is defined by (Aq)~^ = ny{T)cjy-j^ ,^+^-{E), that is, 

A^(r,^^ = 100 MeV) = 2-10^ ( ^^ ^""^ ^ cm. (4.13) 

Using this we may now determine the energy loss. Eg, into the sterile neutrino tower. 
Averaging over the thermal energy distribution of the neutrino, we find 

^ _ 1 mv Jo"" 5n{E,TUE)E^e-^l^^dE 

where the first factor of 1/2 comes because only resonates and < > is the mean 
energy of the neutrinos produced by annihilation. From eq. ([4. 91) we have 



4m2 V / m2 WlOOMeV\ , 



and from eq. (|4^ 

6n = {5p/2p)mnR, (4.16) 
6p/p 

' core ^ \ 

E ^ 



X'{T,E,) ^ 1^ / 30 MeV \ Y J^^ EUE,, + E,fe-^/^dE^ y' 

W 2 \ T J yj;;^ E^e-E^/T dEj, {100 MeV f J 



run = 2(V^)^/^ ~ 10^ eV 



100 MeV 



where the first factor of 1/2 in comes from the fact that the neutrino often elastically 
scatters so that it moves in a direction along which the density does not change. Putting 
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all this together gives 



(100 MeV)^ E^e-^/^ / lOO MeV y^^ 



and hence 



- 10^° ( - ^ e-Me(l/fcT-l/fc 30 MeV)!^ Q u Ig) 



Evaluating the integrals numerically results in 



30 MeV\ 



2.5 



B ^ 0.5 \—^) (4-19) 
so finally 

^5 ~ 4 • 10^^ ( ^^ ^^^ ^ g-Me(l/fcT-l/fc 30 MeV)!!^ _ (4 20) 

Requiring that this be bounded by 3 • 10^^ ergs implies 
m'R _j ( T V ' 



eV 



< a . ^f)-" ( p/^e(l/feT-l/fc30MeV) 



4--2.1 Self-limiting of level-crossing processes 

In fact this bound is too strong because in deriving it we assumed Ui,^ was negligible. If the 
bound is violated then, on a timescale much less than one second, a significant number of 
muon anti-neutrinos oscillate into sterile neutrinos. However the muon neutrinos remain 
as they do not undergo resonant conversion provided is negative. Thus a non-zero value 
of Uy^ builds up driving to zero and thus stopping the resonant conversion!^ One may 
readily determine the energy loss before this happens because the change dn^^j^ needed is 
n„/4 — ny^/2 = 0.13 A'^b- Since the muon neutrinos are pair produced, an equal number 
of muon antineutrinos must have converted into sterile neutrinos. However since the 
antineutrinos are not degenerate and, as discussed above, they have a thermal distribution 
before resonant conversion, they carry approximately 3r energy. Thus the energy loss is 
approximately 1 /7th of the energy stored in the electron and neutrinos sea in the core of 
the supernova, probably just acceptable. 

Does this mean there is no constraint on the muon (and tau) mixing parameter ml 
This depends on whether the mechanism setting to zero is completely efficient. In the 
first half second half of the leptons leave the core [^]. If eq. (t4.21 ) is not satisfied this 



reduction will be made up by a further S'n^^ = buy^j^ Uy^/i: = 0.02 n^, a relatively 
minor correction. More important is the fact that diffusion processes might spoil the local 
cancellation of V^. A crude estimate of this effect may be obtained by noting that the 



A similar effect has been considered in Ref. 
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neutrino cannot be localised within its mean free path, Amfp- As a result, even though 
when averaged over a mean free path the average value, < >, is zero, it may not vanish 
locally. Thus one may expect resonant conversion of to occur locally in regions in which 
is positive and conversion of 17^ in regions in which is negative. Together these will 
not change the lepton number and hence leave < > zero. A naive estimate of the 
sterile neutrino production due to this effect may be obtained if one notes that over the 
mean free path the nucleon density changes by an amount given by 5n„/n„ = Amfp/?'core- 
Hence locally ^ \^Gnn Xmip/'^f core- The overall energy loss is proportional to ^/6V^ 
and hence the rate for resonant conversion when < >= is reduced by the factor 
\/ Amfp / '"core Compared to that given in eq. 

Q. The mean free path is determined by 
the elastic scattering off nucleons via eq. (|4.4|). Using it we find the revised bound 



^ < in-5 r pMe(l/fcT-l/fc 30 MeV) 99^ 

eV ^ ^ VSOMeVy' ^ ^ ' 

This applies to the brane bulk mixing of both muon and tau neutrinos. Notice that if the 
overall neutrino distribution were able to locally compensate the variation of the neutron 
density the bound would become weaker. Analysis of this possibility requires constructing 
a generalised diffusion equation - work on this is in progress. 



4.3 Incoherent production from z/g 

The case of I'e is somewhat different from and Vj- because the pathlength before 
it is eliminated by the process i/g -\- n — > e + p is much shorter than the pathlength for 
i/V annihilation which is relevant to the case of i'^ and Ur- As a result the fe elimination 
pathlength is approximately the same as the elastic scattering pathlength. Hence, there 
is no enhancement factor corresponding to eq. (^]^) and the bound is given by eq. ( [4. 71) 
with r/r = 1, leading to a bound some three orders of magnitude weaker. 



4.4 Coherent production from fg 

The analysis of the case of resonant conversion follows the general lines discussed above 
for the mu and tau neutrinos. In this case 

Ve = V2G{2n^^ +ne + n^^ + n^, - 0.5n„). (4.23) 

An immediate difference is that Ve is positive and somewhat smaller, namely ~ 3 eV. 
This means that it is the that undergoes resonant conversion and the process stops when 
Ve is zero. This corresponds to 5njj^ = —{2ny^ + rig — 0.5n„)/2 = — O.OGn^ corresponding 
to only l/16th of the total energy. Subsequently diffusion causes half of the lepton number 
to leave the core in the first half second and to bring V^ back to zero requires resonant 
conversion of with S'rijj^ = rin/S = — 0.09 ub giving a total loss slightly more than 
l/8th of the total energy. Given that diffusion continues for several seconds the losses 
may become unacceptably large, but a detailed simulation is necessary to resolve this 
question. Assuming that the loss is acceptable one must finally consider the losses due 
to the fact Ve is not locally zero. To determine this we note that the relevant production 
process is e+e~ — > I'eT^e rather than the more rapid process e + p ^ Ve + n. As discussed 
above, conversion occurring with < Ve > zero conserves lepton number and adjusts itself 
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so that the conversion rate for Vf. and are equal. Thus it is sufficient to determine 
the conversion rate for I7g. The energy of the Vf, is only slightly larger than 3T and we 
will take it to be 100 MeV in our estimates. The annihilation rate is dominated by the 
process Ve^V ^ e"*" + n and the corresponding annihilation pathlength is \'(T,Ei;^) = 

30 ^IMM^^ QYn. The resulting energy loss is given by 

dEy 2 2 dE^ ' ^ ^ 

where / '^'^"i^'^^ dE = Qe{T) is given by eq. ( [4.3D , 5n = ^-^^^^ninR, and 7 is taken from 
eq. ( |4.9| ). We estimate m„ for the case < 14 >= in the manner described above giving 

rrin = Y^Amfp/rcore 10^ ( loo^MeV ) Since the neutrino chemical potential is much less 

than that of the electron the energy of the electron and positron will roughly be shared 
between the neutrino and the antineutrino. For the calculation of 6n and 7 we use an 
average neutrino energy of 150 MeV. Using this to evaluate Es leads to the bound 

5 An explicit orbifold model 

In this section, we consider a model with small mixing between SM and bulk neutrinos. 
As a consequence, this mixing can be treated perturbatively and it is unlikely that bulk 
neutrinos play a relevant role in the conversion of atmospheric v^. On the other hand, a 
small bulk component in the electron neutrino might be responsible for the depletion of 
solar neutrinos. Such a possibility, first considered in Ref. [O], has been investigated in 



Ref. 1 22 1 in the context of a model with Dirac bulk masses ||2^. A detailed phenomeno- 
logical analysis has shown that one may go from the limit in which only a single bulk 
state is involved in the resonant conversion to the case where a continuum contributes. 
The former gives the same phenomenology as the standard single sterile neutrino case 
and is disfavoured when the SuperKamiokande (SK) data on the recoil electron spectrum 
is taken into account. However, if more bulk states contribute, the deviation from the 
observed spectrum reduces and a very good fit of the data is possible for a level spacing 



in which three KK modes contribute significantly to the resonant oscillation. In |22] it 
has also been shown that the model can be extended to describe all present indications 
of neutrino masses, including maximal mixing for atmospheric neutrinos in a natural way 
and the LSND signal. In this section, we will show that the same phenomenology can be 
reproduced in the context of a Z2 orbifold model by means of Majorana-type bulk mass 
terms. 



5.1 General structure of the model 

We consider a Z2 orbifold model or equivalently, a model with Z2 invariance. The 
Z2 symmetry is generated by the action of P5 defined in eq. 

(HI)- It 

can be seen from 

Table ||, that only Maj orana bulk masses are allowed by the Z2 symmetry. Also, only 
the brane-bulk mixing terms m, connecting the SM neutrinos to the components of the 
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fermion fields, are allowed while rrf = 0. In addition, we assume that the vector-like 
Majorana term vanishes, that is My = 0. This follows, if the mass Lagrangian respects 
the underlying five-dimensional Lorentz invariance. Alternatively, vector-like Majorana 
masses generated by the VEV of a 5-dimensional vector field are forbidden if one assumes 
that the representation of the gauge group does not mix \I' and Furthermore, for ease 
of notation, we write M = Ms for the Majorana mass. 

Since M is the only bulk mass matrix, it is convenient to consider a basis for the bulk 
fields in which this matrix is diagonal, real and positive, that is M = diag(Mi . . . M^)- 
We continue to call the brane-bulk mixing in this basis m. Specializing the general ex- 
pression ( ^.26 ) to our case, one can then explicitly sum the series to obtain the mass 



matrix 

EmumijMi V ^ 

M2 + n2/i?2 " " ^ TTRmumij cothinRMi) . (5.1) 



m'. 



This equation clearly shows a see-saw type structure, where the "effective" heavy mass 
matrix Mh = diag(M/^ . . . M^) is specified by 

Mf = (7ri?coth(7ri?A//))'^ . (5.2) 

The heavy masses Mj^ involved in the see-saw formula (^]^) is roughly speaking the lighter 
one of the masses Mj and [ttR)^^. In the limit M/ <^ {ttR)^^ we have, in fact, ~ Mj, 
whereas the opposite limit (ttR)^^ <^ Mj leads to ~ (vri?)"^. These results can be 
easily understood. In the first case, only the lowest bulk modes with n = and mass Mg 
effectively contribute while all states with n 7^ are much heavier. Correspondingly, we 
find an "ordinary" see-saw mechanism with the n = bulk modes as heavy neutrinos and 
corresponding heavy masses Mj. In the opposite case, when (vri?)"^ <C M/, the series ( |5.lD 
is dominated by terms satisfying r? jJ^ < M| which are of order —mumij/Mi. Since 
the number of those terms is O {MiR), Mj cancels. It is replaced by R~^ resulting in a 



"geometrical" see-saw mechanism [10|. 



The mass eigenstates follow from diagonalization of the mass Lagrangian, eq. (3.17). 
Actually, as far as neutrino oscillations are concerned, we only need to diagonalize the 
hermitian matrices m'^m' and MnAdn- Using AinAin is clearly more convenient in the 
present case than using Ain- In fact, since M is the only non-zero mass term, A4nA4n is 



diagonal, as can be easily seen from eq. ( |2.14 ). This implies that the fields rj'^ and in 



eqs. ( |3.16a ), although not being mass eigenstates of A4n, are indeed eigenstates of AiliAin- 



Evaluating eq. ( 3.21 ), the electroweak eigenstates as a superposition of mass eigenstates 
are given by 

Here are the light mass eigenstates of mass m^. They are defined by v'^ = UkiVi where 
the unitary matrix U diagonalizes m'^m', that is, diag(mi, m2, ma) = U^m'^m'U . 
For each n > 0, we can define the Z2 eigenstates 
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where ^f, rjj are even and , r]f are odd. By combining modes with opposite n in 
eq. (5.3), we see that the two Z2-odd combinations do not mix with the SM neutrinos. 



This was expected from the general discussion at the end of Section 2. Moreover, only the 
Z2-even combination 



does mix with the SM neutrinos. This means that only one out of the four degenerate 
states ^^j, rj'^j, v'-ni actually couples to the SM neutrinos. The mass M„/ of all 

these states is given by |^ 



Mni = ^M] + ^. (5.6) 
Analogously, in the zero-mode sector the Z2-even field 

N^i ^ ^01 , (5.7) 

with mass Mqj = Mj mixes with the SM neutrinos, whereas tj'q does not. We can, therefore, 
rewrite eq. ( |5.3| ) as follows 

V ^.fc^fc + V^iVo/ + ^/2 V "^/^ N^,, (5.8) 



Eq. (|5.8D expresses the SM flavour eigenstates fg, f^, i^r in terms of the light eigenstates 
/c = 1, 2, 3 and the tower of heavy eigenstates Nni, n > and is, therefore, the starting 
point for studying the oscillation phenomenology in the perturbative regime. The exact 
formulae with matter effects included are explicitly presented in Appendix ^ for the case 
of one family of SM and bulk fields. 

Had we considered a U(l)-symmetric model with Lorentz-invariant Dirac bulk mass 
term only |^^, we would have found similar results. The flavour eigenstates would still 
be given by eq. ( |5.8| ) in terms of light eigenstates Vk and a tower of mass eigenstates Nni, 
where n > 0. However, those states are now expressed in terms of the original states rj^i 
and ^nl as 

p — i<l>nlri T 4- P^'t'nl ri r 

NqI ~ 7/0/ Nnl ^ . (5.9) 



where the angles (pni are defined by the relation m + in/ R = e^^^^y/x^ + n?/R^. The 
masses of the n-th mode is given by eq. ( |5.6| ) but with M/ replaced by the eigenvalues /// 
of the Dirac mass matrix. A crucial difference is that the light states v are massless for an 



^^In Ref. an "effective" Majorana mass term was considered that was induced by a Scherk-Schwarz 
type compactification. In the context of our models with expHcit bulk mass terms such a spectrum can be 
obtained by ve ctor- like Majorana bulk masses and all other mass terms vanishing. This can be explicitly 



seen from eq. (3.14) 



26 



unbroken U(l) symmetry. In a general context, we have already observed this property in 
the previous Section. The mixing matrix in eq. (5.8), therefore, acquires physical meaning 



only once the U(l) symmetry is broken, for example due to small U(l) breaking brane-bulk 
masses or to small bulk Majorana masses. 

Returning to the Z2 model, we should analyze the validity of the perturbative diag- 
onalization that led to eq. ( [5.81 ). Applying the discussion of Section |3.2| , in particular 
eq. ( ^3.16 ), we find that 



£2 = ^ coth(^i?M/) , (5.10) 

should be smaller than one in order to be in the perturbative regime. Then the heavy 
eigenstate components of each SM neutrino flavour i = e, fj,,T in eq. ( |5.8D are small. 
Therefore, the oscillations of SM neutrinos mainly involve the three light Majorana states 
k = 1,2,3 and the corresponding mixing matrix U can be identified with the MNS 
matrix. As discussed above, for such models, some of the observed oscillation phenomena 
can be attributed to oscillations into the three light Majorana neutrinos (playing the role 
of the three standard neutrinos) while some others can be due to oscillation into bulk 
states. 

Since we consider the case that the sterile states play a relevant role in solar oscilla- 
tions only all other oscillation phenomena should be due to oscillation between the light 
Majorana states. We will see that we can naturally obtain Vr oscillation with 

maximal mixing to account for the atmospheric neutrino deficit. At the same time, we 
obtain u^, <-> oscillations with a Am^ larger than the atmospheric one, that can be used 
to accommodate the LSND signal, and two mainly SM neutrinos with masses in the eV 
range, that can provide a cosmologically significant dark matter component. We start with 
an analysis of solar neutrinos in the simple one-family case which we will, subsequently, 
embed in the full three-family model. 

5.2 Solar neutrinos 

We first assume that Ve only mixes with a single bulk fermion ^'g with mass Me through 
a brane-bulk mass term mg. Specializing eq. ( |5.8| ), the electron neutrino can be written 
as the superposition 



where m has been taken positive without loss of generality. The states Nn, where n > 0, 
have mass M„ ~ + n'^/R'^, whereas the mass of the light state is given by 

m'e ^ TTRml coth(7ri?Me) . (5.12) 

Notice that the condition ^ 1 characterizing the perturbative regime is equivalent to 
m'e ^ Me in this case. 



The corrections to eq. ( 5.11 ) and to the mass eigenvalues, as well as the behaviour in the 



presence of matter effects can be obtained by using the exact equations in Appendix 
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For completeness, in Appendix B.l, we have also given the exact equations for the model 



with U(l) lepton number symmetry as briefly discussed in the previous Subsection. From 



eqs. ( [B.6[) we see that the Z2-odd states ^~ for n > 1 and for n > 0, are decoupled as 
expected, a consequence of the exact Z2 parity symmetry. Moreover, from an analysis of 
eq. ( B.22|) , the structure of the mass spectrum in presence of matter effects turns out to 



be the following. 

• There exists one "light" Major ana eigenstate v with mass A' in the range < A' < 
Mf>. In vacuum, A' is approximately given by the see-saw formula ( ^.12 ) implying 
that A' ~ m'e < M^. 

• We have one "zero-mode" eigenstate A^'o with squared mass Aq in the range Mg < 
Ag < Ml + {l/Rf. In vacuum, 

Ao ^ Ml + 2ml - Ml . (5.13) 



For each n > 1 there are two modes with masses A^^ in the range Ml + in/R)"^ < 
Xn± < Ml + {{n + 1) / R)'^ . Their masses in vacuum are given by 

A2±~m2 + — + 2^^^^==4==^~m2 + -— . 5.14 

^2 n ^Ml + n^/R^ ^ R^ 

We see from the last equation that the two n-modes are almost degenerate so that their 
linear combination Nn in eq. ( 5.111) can be considered as approximate mass eigenstates; 



the orthogonal combination is approximately decoupled. The resonance involving the bulk 
state Un becomes relevant for the oscillations of a solar if its energy Ei, is such that 



Ae = 2E^Ve ~ Ml . (5.15) 

Here Vg is the matter induced potential given in terms of the electron and neutron densities 
Ue, rin by Ve = \/2GF{ne — nn/2). The energy-width {5E)n of the n-th resonance, namely 
the width of the energy range where the mixing with the ?2-th bulk states is large, 
can be precisely defined by the condition siv?20n > 1/2. It then follows that {6E)n — 
2y/2me/MnEn, where En is the resonant energy. Notice that, when expressed in terms of 
the parameter y/A^, the width is independent of the resonance, that is {5^/Al)n — \/2me. 
Since {6E)n ^ En+i — En, the resonances take place in a small portion of the energy 
range. This fact, which has already been used in Section ^, considerably simplifies the 
phenomenological analysis. When the neutrino energy is far from the resonant energy, 
the electron neutrino is mainly made of the light mass eigenstate {> if Ag < Ml and it is 
mainly made of the n-th mass eigenstate A„ if 

Ml + ^<A<Ml+^^^. (5.16) 

This determines the behaviour of the electron neutrino survival probability. When the neu- 
trino energy is so small that the value of A^, in the core is smaller than the first resonance, 
that is Ei^ < MI/(2V^°^^), the neutrinos predominantly consist of the light mass eigenstate 
1) and the survival probability is close to 1. For energies larger than Ml/(2V^°^^), on the 



28 




E (MeV) 



10 12 
Erec (MeV) 



Figure 1: Survival probabilities for ~ 2.1 -10 eV and (1/i?, m-e) = (20,0.058) • 
10-3 eV (dotted line), {1/R,me) = (5,0.035) • lO^^ eV (solid line), and {l/R,me) = 
(0.5,0.012) • 10-3 eV (dashed line). 

Figure 2: Recoil energy spectrum associated with the three survival probabilities in Fig. [J. 



other hand, the produced neutrinos are almost orthogonal to the detected ones and the 
survival probability in the adiabatic limit is small. The energy dependence of the solar 
neutrino depletion therefore determines Mg. The Majorana mass Mg must be in fact in 
the range Me = (2 - 3) • 10-^ eV, in such a way that solar neutrinos in the Beryllium energy 
range undergo resonant conversion but pp-neutrinos do not. When the neutrino energy 
increases, the adiabatic approximation fails, the level crossing probabilities grow and the 
survival probability between two subsequent resonances also grows. The rate of growth 
between resonances depends on the mixing in vacuum between the electron neutrinos and 
the mass eigenstates which, for given Mg and R~^ is determined by the brane-bulk mass 
me. On the other hand, the survival probability falls after each resonance, since the num- 
ber of levels to be crossed increases. This behaviour is shown in Fig. 0, where the survival 
probability P(fe — > i^e) is plotted as a function of the neutrino energy for Me = 2.1-10^3 gv 
and three values of 1/R and nie. The three values of l/R, 0.02 eV, 0.005 eV and 0.0005 eV, 
represent a sample of the different possibilities offered by the model. The values of the 
brane-bulk couplings have been determined in each case by fitting the total rates measured 
in the SuperKamiokande [0, Homestake |42| and Gallium |43, 0, 45| experiments. 
They are nie = 0.58 • 10-^ eV, rUe = 0.35 • 10-^ eV, = 0.12 • 10-^ eV respectively. The 
corresponding predictions for the recoil energy spectrum in SuperKamiokande are shown 
in Fig. ^ together with the latest data |2£]. 
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The dotted line corresponds to the regime M^R <C 1. The probabihty does not depend 
on the precise value of 1/R for the energy range shown in Fig. |l] as long as 7Me. In 

this case, in fact, 2E^V^°'^^ < + R? , so that the resonant mixing between the electron 
neutrino and the sterile tower is significant only for the lowest state. The mixing angle 
with the lowest state Nq is — '^el^'^e- The value of mg corresponds to sin^ 20o — 3- 10~^ 
and to a light neutrino mass of 1.7 • 10~^ eV. The phenomenology of this case is the same 
as the one for a model with a single sterile neutrino — the Kaluza Klein origin makes 
no difference. The SuperKamiokande collaboration has claimed that their data on the 
day-night asymmetry and especially on the electron recoil energy spectrum disfavour such 
a scenario involving oscillations into a single sterile neutrino. This is confirmed by Fig. 
where the dotted line has the highest slope and therefore gives the worse fit of the data0, 
which is well compatible with a flat spectrum. 

The presence of additional sterile states reduces the slope of the predicted energy 
spectrum. This is illustrated by the dashed line in Fig. |l] and the corresponding dashed 
histogram in Fig. |2|. They both correspond to the small XjR case = 0.5 • 10-3 eV m 
which Mgi? > 1. This situation is quite new because a relatively large number of levels 
undergo resonant conversion. The single resonances are not visible in the figure due to the 
average over the neutrino production point in the sun. The mass of the light neutrino is 
now 0.9 • 10"^ eV. This case leads to a significant improvement with respect to the single 
sterile neutrino case with the fit of the total rates (the energy spectrum) at 65% CL (50% 
CL). Finally, the solid lines in Fig. Fig. ^ represent the intermediate case 1 jR = 0.005 eV 
in which, besides the light state three bulk states, Nq, Ni and participate in solar 
oscillations. Fig. |2| shows that the values of the bulk and brane-bulk masses giving the 
best fit of total rates (within a 25% CL) also beautifully fit the energy spectrum (10% 
CL). Besides the three examples discussed, the whole range of possible values oil/R also 



includes the situation in which M^R ~ 1 [22|, that approximately reproduces the case 



considered by Dvali and Smirnov |12]. 



Unlike the energy spectrum, the day-night asymmetry ^] does not discriminate 
between the different possibilities discussed. In fact, the day-night asymmetry turns out to 
be small compared to the experimental errors in all cases considered above and therefore 
is in reasonably good agreement with the present data, given that SK finds only a 1.2 o" 
deviation from zero asymmetry. 

For our three sample cases the values of m^R/ eV are of the order of 10^'^ and, hence, 
compatible with the supernova bounds derived in Section ^. 

The scenario discussed leads to various experimental signatures. First of all, a crucial 
test will be provided by the neutral/charged current ratio measurement in SNO [47|. 



Moreover, different sterile neutrino scenarios can be distinguished by the shape of the 
energy spectrum. As discussed, the present data already favours a flatter spectrum, which 
can be obtained if more than one sterile neutrino takes part in the oscillations. The charged 
current spectrum that will be measured by SNO could provide additional information. 
Finally, the measurement of a day-night asymmetry signiflcantly lower than zero would 
disfavour the scenario with many sterile neutrinos. 



^^With free spectrum normalization, such a fit is unacceptable at 95%. A better agreement can however 
be found for mixing angles at the lower border of the range allowed at 90% CL by the total rate fit |22|]. 



30 



5.3 Atmospheric neutrinos 



Let us now discuss atmospheric neutrinos. As already pointed out, in a model with 
small (perturbative) brane-bulk mixing, the small bulk components of the SM neutrinos do 
not contribute significantly to atmospheric neutrino oscillations. Hence, these oscillations 
must be due to the three light mass eigenstates. We can, therefore, focus on the first part 



of eq. (5^) given by 

Ui= ^ Uikh ■ (5.17) 



k=l,2,3 



The light masses and the MNS mixing matrix U are determined by the light mass matrix 
m' in eq. ( |5.1D . Before constructing an explicit model, let us first discuss some phenomeno- 
logical requirements on this mass matrix suggested by atmospheric neutrino oscillations. 
The simplest texture leading to maximal Vn <-> i't mixing is 



m' = I 1 ) m^r , (5.18) 

r 1 

where the e parameters are much smaller than one. The parameter a is not necessarily 
small but it will be in our explicit model described below. 

In the limit in which the e parameters are vanishing, the muon and tau neutrinos are 
superpositions of two degenerate states 1)2 and 0^ with mass |?7t-^t-|. The small e parameters 
(assumed to be real for simplicity) are necessary to generate a mass splitting 

Am23~e|m^^|, e = + (5.19) 

between the degenerate states and therefore an "atmospheric" squared mass difference 
Am^rpjyj = Am|3 ~ 2e\m^r\'^. The corresponding mixing angle 6^r is almost maximal and 
it is given by sin^ 20^,- ~ 1 — (e^ — er)"^ /A. For Am^g — ^"^13 ~ ^"^atm parameters 
Ce/i; Cer are constrained to be small by the i^e disappearance experiments. 

Notice that the simple mechanism described to generate a maximal z/^ <-> i^r mixing 
does not work in a three neutrino scenario aiming at explaining at the same time the solar 
neutrino data. This is because the other two squared mass differences available turn out 
to be larger than the atmospheric one, unless the parameter a in ( ^.18 ) is close to 1 and 



all neutrino masses are degenerate [^]. If this is not the case and all e < 1/2, we have, in 
fact, 

Am^i ~ ^"1-31 ^ I^T-itrP > 2e|m^^p = Am^rpj^ . (5.20) 

As a consequence, there is no room for the small Am^ required by the solar data. In our 
case, of course, such a small Am^ is not needed since the solar neutrino problem is solved 
by oscillations into bulk neutrinos. 

If the texture ( 5.18| ) accounts for atmospheric neutrino oscillations, the two degenerate 



neutrinos can provide a cosmologically significant source of dark matter. Moreover, the 
electron neutrino can oscillate into muon or tau neutrinos with a squared mass difference 
Am'^ > Am^rpj^ and small amplitudes 4eg^, Ae^^ respectively. If the parameter a is small, 
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we have Am ~ m^r- It is obviously tempting to associate such oscillations with the signal 



observed by LSND. This then determines the size of m^,- to be |m^T-| — y ^"^lsnd- '^^^ 
MiniBooNE experiment will test this possibility and will cover a relevant portion of the 
parameter space for such short-baseline oscillations. A short-baseline neutrino factory 
could further extend the sensitivity in the mixing parameter 



Let us now see how the texture ( 5.1^ ) can be obtained in a model in which the sin- 



gle family discussion of solar neutrino oscillations given above can be embedded. We 
work in the context of the Z2 invariant models with Majorana bulk mass described in 



Subsection 5.1 



A simple way to obtain the zeroth order form of eq. ( |5.18| ) in which a and the e 
parameters vanish is to assume three bulk fermions ^^i, ^r, a diagonal brane-bulk 
coupling m and a Majorana mass matrix M with two off-diagonal non vanishing entries. 
Concretely, 



m 





































/ 



(5.21) 



where M^,- can be made real by a phase-redefinition of the bulk fields We note 

that, for a direct application of the results of Subsection |5.1| to this case, we should first 
diagonalize the matrix M and rewrite the brane-bulk coupling m in the corresponding 
basis. Alternatively, one can use eq. (P^) with m and M as m eq. (|5^ . Then, we find 



from eq. ( ^.1| ) that the light mass matrix m' has indeed the structure ( ^.18 ) with vanishing 
parameters o, e and the overall scale m^r given by 

m^r ~ —TTRm^mr coth^TrRM^r) ■ (5.22) 

The bulk components of and Vr are small if 

TTR\m^\'^ coth{7r KM f,r) <^ M^r and 7ri2| ^ coth(7r i?M^^) <C M^^ , (5.23) 

which implies that 

\m^r\<^M^r- (5.24) 



The form ( ^.21 ) of the brane-bulk mixing and of the Majorana mass matrix can be 



enforced by means of symmetries acting both on brane and bulk fields. For example, let us 
introduce three individual lepton numbers L^^^^r for each generation with the bulk fermions 
^e,fi,T carrying the same charge as the corresponding neutrinos on the brane. Let us further 
require that — and Lf, are both conserved. This leads exactly to the structure given 
in eq. ( |5.2l| ). In order to generate the small e corrections to the texture ( ^.18 ) one has to 



break these symmetries. This can be done by small off-diagonal brane-bulk couplings or 
small Majorana masses replacing the zero entries in the matrices ( ^.21 ). 



To see this in more detail, let us first discuss how solar neutrino oscillations along the 
lines of Section 5^ can be embedded into the model. As a consequence of eqs. ( ^.21 ), the 



electron neutrino only mixes with the bulk field ^e- The solar neutrino phenomenology 
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described in Section follows once Lg is broken by a non-vanishing Majorana mass 
Me. The smallness of Mg ~ (2~3) • 10"'^ eV compared to M^r ^ l^i^rl > y'^^^'^iTM 
then a consequence of the approximate conservation of Lg . Analogously, a small breaking 
of — Lr can produce the small diagonal entries and Cr necessary to generate the 
atmospheric Am^. This can be accomplished by switching on the — breaking and 
Le conserving entries in the bulk Majorana matrix or in the brane-bulk mixing matrix. 
Once Lg and — are both broken, the v^-Vp^ mixing in ( ^.18 ) can be easily generated 



in the same way. In general, — L-y can be broken both by bulk Majorana mass terms 
and by brane-bulk couplings. The most general case, corresponding to replacing all zero 



entries in the texture (5.21) by small corrections, leads to a quite complicated mass matrix 
for the light neutrinos. As an illustration, let us consider a somewhat simplified case where 
we break — L^- by brane-bulk terms only. That is, instead of eq. ( ^.21 ) we now consider 
matrices 



(5.25) 



where Sme^, dmer and 6m are small corrections. Expanding to first order in these 
corrections and using the parameterization ( |5.18| ), we find for the e perturbations to the 
mass matrix 





( rUe 


SrUep 






m = 




nip 


5m pr 


M 




^ SrUre 


Snirp 


rrir ) 





Me 














Mpr 





Mpr 






26m 



me coth(7ri?Me)(5mg^ -|- m^ cot]i{-K RMpr)6mre 
m^mr coth(7ri?M^T-) 

me coth.(TT RMe)SmeT + rrir coth(7ri?M^,-)5m^g 
m^mr coth.{'KRMpr) 

The atmospheric square mass difference is then related to the corrections by 



(5.26) 
(5.27) 
(5.28) 
(5.29) 



2|m^^|2 V rn„, m-r 



+ — ^ . (5.30) 



The amount of — Lj- breaking needed for a correct account of atmospheric neutrino 
oscillations is, therefore, variable since it depends on m^r which is not fixed by solar or 
atmospheric data. If, in addition, we require an explanation of the LSND signal we have 



^"^LSND resulting in a — breaking e of the order y Am^rpj^/AmLSND- 
For a hierarchical structure of the brane-bulk mass matrix, that is rrig <C m^ <C m^ and 
^im-re ^ 6m pe, the e-/i mixing angle 9ep is given by 



e/i ~l~ (^er ^m^e 



V2 V2'. 



m 



(5.31) 



This angle should then be in the LSND range which can be easily accommodated by 
choosing 6mep- Notice that such a small t'g-t'^ mixing can also generate a small mixing 
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between the electron neutrino and the towers of sterile neutrinos associated with ^^,t, 
heavier than the one we have used to deplete the solar neutrino flux. Due to the large mass 
splitting between the lightest electron neutrino and the lowest mode of the heavy towers, 
such a small mixing can not be MSW-enhanced and consequently does not affect the solar 
neutrino phenomenology. However, it may give rise to a non-negligible contribution to the 
light electron neutrino mass. 

It is clear, that a parallel discussion can be carried out for — Lj- breaking induced 
by bulk Majorana mass terms. However, the above example is sufficient to show that all 
oscillation data can be explained in terms of our model. We, therefore, will not carry this 
out explicitly. 

We have already mentioned that parameters in the model associated with the first 
generation are compatible with the supernova bound. However, we also have to worry 
about the supernova bound due to resonant conversion of and into their respective 
Kaluza-Klein states. As the above equations show, the particle phenomenology of the 
model is not very sensitive to M^^ as long as this mass is larger than m^r and smaller 
than the fundamental scale. Therefore, we can choose M^^ larger than V2EV where V 
is the matter potential in the core of the supernova and E is the neutrino energy. This 
avoids any problem with the supernova bound due to conversion of and i^r into their 
associated bulk states. As discussed, the large value of M^,- (as compared to Mg) can be 
attributed to the approximate conservation of L^. 

To summarize, in the framework of a Z2-symmetric theory we have accommodated solar 
neutrino oscillations, naturally maximal <-> I'r oscillations, two degenerate neutrinos 
providing a cosmologically significant source of dark matter and I'e ^ oscillations with a 
Am^ larger than the atmospheric one, that can be used to accommodate the LSND signal. 
The spectrum of parameters in the bulk Majorana mass matrix and the brane-bulk mixing 



matrix (5.25) we need is the following. 



The value of should be compatible with its lower experimental limit emerging 
from Cavendish-type experiments. 

To get the best agreement with the SuperKamiokande energy spectrum, 1/R should 
be less that 10~^ eV, in such a way that more than a sterile state is involved. 

To explain the solar neutrino deficit has to be broken by an electron bulk Ma- 
jorana mass Me ~ (2 - 3) • 10-3 eV. Furthermore, we need the electron brane-bulk 
coupling to be of the order me = O (iQ-^^"^-' eV). 



If the LSND signal has to be explained the parameter \m^r\ as given in eq. ( 5.22| ) 



should have the value \m^r\ = a/^'^lsnd ~ ^ (l^V). This requires that the quan- 



tity yjm^mr is close to the geometric average of y Am^gND and the smallest between 
M^r and R-^. 

• We need M^-r 3> | m^r | , in order to have a small bulk component in and Uj- and a 
quasi- unitary light mixing matrix U . A value of M^^ larger than ^/2EV avoids any 
potential problems with the conversion of u^^ and u-r into their associated bulk states 
in the supernova. Also the hierarchy between and rrir must be smaller than the 
hierarchy between \m^-r\ and M^r- 
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• To explain the atmospheric squared mass difference we need to break L^ — L-r. If this 
is done be brane-bulk mixing and we use the LSND result the level of this breaking 
is specified by the condition dm^r/m-n — ^"^atm/'^"^lsnd • 

• In order to get the e-/x mixing required by LSND we need to consider terms breaking 
Le and — Lr- If this is achieved by perturbing the brane-bulk mixing the level 
of this breaking is determined by Oe/i = Sm^e/mn,^- A smaller value of 9e^ than 
needed to account for the LSND signal could be within the reach of MiniBooNE or 
a short-baseline neutrino factory. 

Finally, let us mention that there are alternative models within the general class con- 
sidered in this paper which provide similar possibilities to explain the various oscillation 
phenomena. For example, in the context of the models with U{\) lepton number symmetry 
and Dirac bulk masses which we have briefly discussed, an analogous (and complementary) 
scheme using the brane-bulk mixings as a source of f/(l)-breaking and light neutrino 
masses can be constructed Then, the overall U{\) can be considered as a subgroup of 
the individual lepton number groups and the model with — and Lg symmetry arises 
at an intermediate stage in the necessary breaking of these individual lepton numbers. 

6 Conclusion 

In this paper, we have systematically studied five-dimensional brane-world models of neu- 
trino physics based on flat compactifications. Motivated by the expectation from string 
theory, we have particularly focused on the effect of bulk mass terms for the sterile neutri- 
nos. For such models we have presented a number of general results. We have pointed out 
that there are four different types of bulk mass terms, namely scalar Dirac and Majorana 
masses as well as vector-like Dirac and Majorana masses. While the scalar masses can be 
thought of as being generated by the VEV of a bulk scalar field and respect five-dimensional 
Lorentz invariance the vector masses are due to the VEV of the fifth component of a bulk 
vector field and break five-dimensional Lorentz invariance. We have analyzed the (unper- 
turbed) spectrum of bulk Kaluza-Klein states in the presence of these mass terms and have 
found a crucial difference between scalar and vector- like masses. Whereas scalar masses 
constitute a lower bound on the spectrum this is not the case for vector-like masses. For 
arbitrary values of these masses (and in the absence of scalar masses) there will always 
be states with masses of order 1/R. In the presence of small brane-bulk couplings, we 
have presented perturbative results based on a see-saw type approximation. In particu- 
lar, we have shown in general that, in this approximation, the electroweak eigenstates are 
mainly given by a superposition of three light mass eigenstates with small admixtures of 
massive Kaluza-Klein modes. We have also shown that the relevant parameters associated 
with these two sectors can be controlled quite independently by choosing the various mass 
terms in the action appropriately. This opens up the possibility of associating some of the 
observed oscillation phenomena with "standard oscillations" into the light states and some 
others with oscillations into heavy bulk states. To complement the perturbative approach, 
we have also presented a general formalism and results for the exact diagonalization of the 
mass matrix. 

These general results have been applied to construct a specific brane-world model 
compatible with all presently observed oscillation phenomena. The model is based on 
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a Z2 orbifolding and broken Lg and — Lj- symmetries. The solar neutrino deficit is 
explained by small mixing angle oscillations into the tower of heavy bulk states. It has 
been shown that this, contrary to the case of small mixing angle oscillations into a single 
sterile neutrino, is compatible with recent SuperKamiokande results. The model also 
naturally leads to a maximal •s-^ Ut mixing that accounts for the atmospheric neutrino 
result. In addition, it allows to accommodate the LSND observation and may lead to a 
significant neutrino dark matter component. 

We have also analyzed constraints from supernova energy loss on such brane-world 
models and have demonstrated that our specific model is consistent with these constraints. 

In summary, we have shown that the inclusion of bulk mass terms leads to a rich 
structure of brane-world models for neutrino masses which can be used to construct models 
compatible with all known experimental results. 
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A Spinor conventions in four and five dimensions 

In this Appendix, we collect the properties of five-dimensional gammma matrices and 
spinor s which we are using in the main part of the paper. 

The five-dimensional gamma matrices 7^, where a,/?, ••• = 0,1,2,3,4, satisfy the 
Clifford algebra 

{7a,7^} = 2r?a/3 = 2diag(l,-l,-l,-l,-l) . (A.l) 

Note, that we are using the "mostly minus" Minkowski metric throughout. As usual, the 
generators T,aj3 of the Lorentz group <S'0(4, 1) can be obtained by 

^a0 = \ba,lp]- (A.2) 

The hermitian conjugate of the gamma matrices is generated by 70 according to 

(7a)^ = 7o7a7o • (A.3) 

Similarly, to express the complex conjugate and the transpose of the gamma matrices we 
introduce a matrix B and a charge conjugation matrix C. Their defining properties are 

(7a)* = B^^B-^ (A.4) 

{laf = C^aC-^ . (A.5) 

Obviously, this implies that C = B^jq. We note, that the matrices B and C square to 
— 1 unlike their four-dimensional counterparts which square to one. 
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A five-dimensional Dirac spinor \I/ has four complex components. Under infinitesimal 
Lorentz transformations it transforms as 



5* = ^LOap^ap"^ ■ (A.6) 

We define the conjugated spinor ^ and the charge conjugated spinor by 

^ = ^rt^o , = B'^^* = C*^ . (A.7) 

Note that, under infinitesimal Lorentz transformations, transforms in the same way as 
Furthermore, we note that {"^'^y = — The minus sign in this relation is due to the 
above mentioned fact that = = — 1 and precludes the existence of five-dimensional 

Major ana spinor s. 

For our purpose of dimensional reduction to four dimensions, we need to make contact 
between the above five-dimensional spinors and Weyl spinors. To do this, it is useful to 
introduce explicit representations for the gamma matrices. We choose 

- ^, - 



7m = [^^ oj , 74 = 70717273 = ^ J (A.8) 

where = (I2, — Ci), = {I2, o'i) and (Tj are the Pauli matrices. Furthermore, we have 
the matrix 75 defined as usual by 

75 = m = f ° ) (A.9) 



.0 -I2, 

and satisfying (75)^ = 1. For the matrices B and C we can choose, in this representation 

S = 7274 = - (° , C = 727074 = (l (A.IO) 

where e = ia2 is the two-dimensional epsilon symbol. 

A five-dimensional Dirac spinor is decomposed into two left-handed Weyl spinors ^ 
and r] as 

* = ffl • (A-11) 



'7 

The conjugated spinor ^ is defined by ^ = e^*. Using the above explicit representation, it 
is easy to obtain the relations 

* = (r?t -ee) , = , = (-^t _^T,) (A.12) 

which will be useful for the dimensional reduction to four dimensions. To write our four- 
dimensional effective action in a compact form we furthermore adopt the following notation 
for Weyl mass terms 

66 = e6 . (A.13) 

where 6 and 6 are two (left-handed) Weyl spinors. 
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B Exact diagonalization 



We first consider tlie general case of three standard model neutrinos = (z/e,f^,i^- 
and K "bulk" neutrinos x = (xi • • • Xk)'^ ■ Tlie associated Majorana mass matrix is of the 
general form 

T T 
X 

where v is a k x 3 matrix and D is a k x k symmetric matrix. When the set of x fields 
is allowed do be infinite, the results can be applied to the diagonalization of the mass 



matrix (2.18), originating from the five-dimensional brane-world theory, by making the 



obvious identifications 

U = {Ue, Ufj,, Vrf X = Hn, rin)nez^ (B.2a) 
= (m^,m''^)„gZ B = {MnS-m,n)n,m€Z- (B.2b) 

While, in this paper, we focus on cases where D has such a block structure, the general 
formulae below also apply to more complicated cases. A more complicated D is generated 
if the bulk mass terms effectively depend on the coordinates of the additional dimension. 
Such an explicit coordinate dependence may arise in orbifold models or due to a non-flat 
vacuum solution. 

Since, with the usual approximations, the equations for neutrino propagation only 
involve A4^M, we actually diagonalize the matrix 

Ml.^M^M,{^^ ^,%,). ,B.3) 

where the hermitian matrix A = y/2EyV takes matter effects into account. As usual, 
Ey represents the neutrino energy and V the matter induced potential. The eigenvalue 
equations for the generic mass eigenstate 



+ P^X (B.4) 



(in obvious matrix notations) are 

f (v+v + + v"^D/5 = X^a 
[ Dva + DD/3 + = \^ (3 

where is the corresponding real (possibly negative) eigenvalue. 

As shown in Appendix |B.2| and further explained in the main text, it often happens 
that some of the eigenvectors of the matrix DD are also eigenvectors of the matrix M 
and as such decouple from the light neutrinos u. The vector space generated by such 
"unperturbed" mass eigenstates coincides with the maximal space of vectors in the form 



eq. (B.4) such that a = 0, /? is orthogonal to both v and Dv and the space is invariant 
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under DD. In other words the unperturbed eigenstates are in the form i>unp = P^X^ with 
(3 solving the following set of equations 



' DD/? = X^(3 

v^D/3 = 
. = 



(B.6) 



Let us now derive the eigenvalue equations for the "perturbed" eigenstates. In this case, 
it is reasonable to assume that is not an eigenvalue of DD, so that — DD can be 
inverted |^. Then, using eqs. ( |B.5| ) v^D/? and /3 can be expressed in terms of q = 
and a. In particular. 



1 



A2 - DD 



[Dva + v(q;) 



(B.7) 



By taking appropriate projections, one then gets the following equations for a = w'^ (3 and 
a = (A2-vtv-^)-ivtD/3: 



1 



V a 



1 



A2 - DD 

A2 _ vtv - A - v^D 



A2 - DD 



Dva 



1 



A2 - DD 



Dv I a = v^D 



1 



(Bi 



A2 - DD 



va 



Eqs. ( |B.g| ) can have two different sets of solutions. In some cases there can be "decoupled 
solutions" corresponding to eigenstates decoupled from the SM eigenstates, that is a = 0. 
The decoupled states of squared mass A^ are therefore in the form 



~t T 

Udec = a'v 



1 



A2 - DD 



X , 



(B.9) 



where a = v'^/J must satisfy 



1 

v^D 



1 



A2 - DD 

1 



v a = 



(B.IO) 



-va = 



A2 - DD 

Although such decoupled (but not unperturbed) eigenstates are not always present, we 
will see in Appendix B.l an example in which they are relevant. Also note that the 
equations for both the unperturbed and the decoupled states do not depend on matter 
effects encoded in A. Therefore, these states are insensitive to matter effects as one would 
intuitively expect. 

The generic solutions of eqs. ( p.8| ) correspond to states which couple to the SM eigen- 
states u and, therefore, have a ^ 0. In this case, if follows from eqs. ([B.8|) that A^ is an 



'^^The eigenvalue of a perturbed eigenstate could "accidentally" be an eigenvalue of DD as well. Here, 
we do not consider such pathological cases explicitly. 
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eigenvalue if the following equation has a non-trivial solution for a: 



1 - V 



1 



A2 -DD 



a 



V 1 - 



1 



1 



-Dva . (B.ll) 



A2-DD V A2-DD y A2-DD 

In the case of a single family and no matter effects this eigenvalue equation becomes 



A^ 



1 



A2 - DD 



v^D 



1 



A2 - DD 



(B.12) 



The eigenspace corresponding to the eigenvalue A^ is parameterized by the values of a 
solving the above equation. By using eq. ( |B.7| ) with a obtained from the first of eqs. (B^ 
we find the following mass eigenstate associated with A^ and q. 



+ v^D 



1 + 



1 



A2 - DD 



V 1 - 



1 



A2 - DD 



-1 



1 



A2 - DD 



X 



(B.13) 



We now specialize these results to the case of a mass matrix D with block structure 
as in eqs. (|B.2D , a single SM family v and a single bulk fermion (^n)^n)nez- This is done 
for two cases, namely for a Dirac and a Majorana bulk mass term. 

B.l Dirac bulk masses 

We first consider the simple case with vanishing Majorana-type terms, that is Ms = My = 
and m'^ = 0, and a Dirac bulk mass term with a scalar components only, that is // = //g, 
/iy = 0. Such a pattern is realized in the context of a model with U{1) lepton number 
symmetry. Furthermore, we also include matter effects. 

The unperturbed eigenstates can be derived from eqs. ( |B.6D . For each n > 1 there are 
two unperturbed eigenstates with mass y^/i^ + jB? ^ namely 



and 



V2 



(B.14) 



Here, the phase is defined by ji + in/R = e**^" iJ? + jB? . 

In this model decoupled eigenstates are present as well. In fact, as a consequence of 



the U(l) symmetry, the second equation in (B.IO) is automatically satisfied. This is a very 
particular property of the ?7(1) symmetric model and will generically not be the case for 



other model. The first equation in (B.IO) then leads to the eigenvalue equation 



7ri?|m|^ cot ttR^/ A^ — 



For each A solving this equation one obtains a decoupled state of the form 



-y 

/at 



in ; 



(B.15) 



(B.16) 
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The normalization factor in eq. ( B.16| ) is given by 



2N 



A2 



+ 



A2 



(l + 7r2i?2|m|2) 



|m|' A^ — 

The equation for the masses A of the coupled eigenstates is given by 



(B.17) 



A^ 



: cot T^RyJX^ - /u2 + A = \^. 



(B.18) 



This eigenvalue equation has a single solution with A^ < n^, as can be seen by using the 
relation coi{'\/ —x) / \J —x = — coth(y^)/y^. In vacuum, we have A^ = 0, as a consequence 
of the U(l) symmetry. The associated state with mass A^ reads 



ne2 



m{pL — in/R) 



(B.19) 



The normalization factor N in eq. ( [B.19| ) is given by 



2N 



{y-AY + 2A\my \\l + TT^R^\m 



X^\m\ 



+ 



A 



A2-^2 



(B.20) 



In vacuum, the masses of the coupled eigenvalues coincide with the masses of the decoupled 
ones, with the exception of A = which is only a coupled eigenvalue. Also, the normal- 
ization factors in eqs. ( pl^ ) and (pl9D coincide in vacuum. Coupled and non-coupled 
states can be combined to form left and right-handed components of Dirac fermions. 



B.2 Majorana bulk masses 

We now consider the case of a vanishing Dirac mass, = and furthermore assume that 
mP = and My = 0. We keep the scalar Majorana mass M = Ms- This theory is 
invariant under the Z2 symmetry. 

The unperturbed eigenstates can be derived from eqs. ([B.6|). For each n > 1 there 



are two unperturbed eigenstates with mass + jR? corresponding to the Z2-odd 

combinations 

r]Q is the unperturbed eigenstate corresponding to n = 0. There are no decoupled eigen- 
states. 

The exact masses A of the coupled states are determined by eq. ( |B.11| ), which in our 
case becomes 

[X'^{l-x)-A]{l-x) = M^x^, (B.22) 



In vacuum the normalization of the light state is A'' = 1 + 7r7?|mp coth{TrRfi)/ fj.. 
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where 



X = TrRlml 



,cot 7ri2VA2 - M2 



VA2 - M2 



By solving for x, Eq. (B.22) can be written as 



cot ttRVX^ - M2 



A 



VA2 - M2 X±M 
The corresponding mass eigenstates are 

1 



E 



nG2 



A2-M2-n2/i?2 V^R ' 1-x 



.n M 



where is a normahzation factor. 



(B.23) 



(B.24) 



(B.25) 
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